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Abstract 

This  research  presents  a  method  using  modal  analysis  by  which  electromagnetic 
characterization  of  materials  in  a  partially  filled  rectangular  waveguide,  having  a  sin¬ 
gle  top  air  gap,  can  be  accurately  performed.  Thermal  expansion  of  waveguides  com¬ 
monly  occurs  during  high-temperature  measurements,  resulting  in  air  gaps  between 
the  sample  and  waveguide  walls.  Higher  order  modes  are  excited  by  the  discontinuous 
geometry,  which  are  not  accounted  for  in  most  closed  form  extraction  algorithms.  A 
correction  must  be  applied  that  considers  the  complex  power  transmitted  and  stored 
by  higher-order  modes,  not  merely  the  dominant  mode.  Characterization  independent 
of  sample  distance  from  the  calibration  plane  is  also  presented. 

Expanding  upon  previous  analysis  of  partially  filled  rectangular  waveguides,  a 
modal  solution  for  a  single  air  gap  between  the  top  of  the  material  sample  and  the 
waveguide  wall  is  developed.  The  analysis  is  performed  on  samples  of  dielectric  to 
verify  the  method,  and  further  tests  are  performed  on  magnetic  shielding  material. 
Boundary  conditions  between  the  empty  and  partially  filled  regions  are  formulated 
so  it  is  only  necessary  to  explicitly  satisfy  two  of  the  existing  three  field  components, 
the  third  being  linearly  dependent. 

Calculation  of  the  complex  permittivity  and  permeability  of  magnetic  shielding 
material,  within  10%  of  the  true  value,  was  achieved  by  using  less  than  20  modes. 
Measurement  improvement  of  real  permittivity  is  the  strongest  feature  of  the  algo¬ 
rithm.  Modal  corrections  of  samples  converged  slowly  or  not  at  all  if  a  large  air  gap 
(>  90mils)  was  present,  indicating  the  sensitivity  of  the  solution.  The  inclusion  of  the 
calibration  plane  independence  analysis  greatly  improved  the  level  of  performance  of 
the  modal  correction.  Results  are  presented  in  S-band  and  X-band,  although  the 
concept  is  applicable  to  all  frequency  bands. 
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MATERIAL  CHARACTERIZATION  IMPROVEMENT 


IN 

HIGH  TEMPERATURE  RECTANGULAR  WAVEGUIDE 

MEASUREMENTS 

I.  Introduction 

Electromagnetic  characterization  of  materials  quantifies  the  response  to  applied 
electric  and  magnetic  fields  through  the  use  of  the  complex  quantities  permittivity  (e) 
and  permeability  (/i),  respectively.  The  real  part  of  each  complex  number  is  related 
to  energy  storage  in  the  material,  while  the  imaginary  part  accounts  for  conduction 
losses,  manifested  as  thermal  energy. 

In  applied  electromagnetics,  the  parameters  e  and  p  are  routinely  sought  by 
microwave  engineers,  since  these  parameters  determine  the  behavior  of  fields.  Mate¬ 
rials  can  be  either  non-magnetic  dielectrics  or  magnetic  materials  with  both  electric 
and  magnetic  losses.  It  is  common  for  certain  types  of  materials  to  be  used  to  ab¬ 
sorb  and  dissipate  electromagnetic  field  energy,  such  as  that  transmitted  by  a  radar. 
These  shielding  materials  are  effective  at  reducing  the  echo  area  or  radar  cross  section 
(RCS)  [20].  The  materials  used  to  achieve  this  reduction  typically  vary  in  weight, 
composition,  effective  frequency  band,  and  durability,  as  well  as  method  of  applica¬ 
tion. 

A  highly  relevant  topic  in  RCS  engineering  is  the  reduction  of  scattering  from 
the  exhaust  cavity  of  aircraft  engines.  Shaping,  the  primary  tool  for  RCS  reduction, 
cannot  be  used  on  existing  legacy  aircraft.  Therefore,  radar  absorbing  shielding  ma¬ 
terial  must  be  applied  to  needed  areas.  It  is  necessary  that  the  absorber  be  able 
to  withstand  the  extremely  high  temperatures  of  the  exhaust  environment,  typically 
2500°  F,  without  performance  degradation.  Ceramic  shielding  material  is  commonly 
used,  because  of  its  durability  under  heat. 
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As  mentioned  before,  knowledge  of  e  and  //  of  the  shielding  material  is  necessary 
to  determine  the  response  of  the  electromagnetic  fields  and  the  effective  RCS,  that 
is,  whether  the  shielding  material  will  meet  the  requirement  for  RCS  reduction.  Mi¬ 
crowave  measurements  are  used  to  calculate  the  electromagnetic  parameters  of  heat- 
resistant  ceramic  shielding  material,  but  must  be  performed  at  elevated  temperatures 
to  properly  simulate  the  performance  environment.  This  requirement  introduces  sev¬ 
eral  complications  into  the  measurement  process,  which  would  otherwise  be  trivial. 
Thermal  expansion  of  the  metal  waveguide,  which  is  normally  greater  than  that  of  the 
ceramic  shielding  material,  introduces  air  gaps  between  the  sample  and  the  waveguide 
walls,  leading  to  the  excitation  of  higher  order  modes;  this  situation  will  be  referred 
to  as  a  partially  filled  waveguide  (PFW).  If  a  gap  is  nonexistent,  then  the  situation 
is  a  fully  filled  waveguide  (FFW). 

The  use  of  closed  form  algorithms,  such  as  Nicolson-Ross-Weir  (NRW),  does 
not  account  for  the  complex  power  loss  to  higher-order  and  evanescent  modes  [23,27]. 
These  must  be  considered  if  an  accurate  extraction  is  desired.  Occasionally,  due  to 
the  extreme  heat  or  imprecise  placement,  the  sample  will  shift  longitudinally  away 
from  (or  towards)  the  calibration  plane,  introducing  a  phase  shift.  If  the  phase  shift 
is  not  compensated,  an  extraction  of  the  electromagnetic  parameters  will  be  flawed. 
This  thesis  uses  modal  analysis  techniques  to  overcome  the  problem  of  an  air  gap  and 
is  combined  with  a  new  method  for  reference  plane  independent  measurements  in  a 
rectangular  waveguide  containing  magnetic  material. 

1.1  Problem  Statement 

A  need  exists  to  accurately  characterize  magnetic  shielding  materials  at  ex¬ 
tremely  high  temperatures.  During  high  temperature  tests,  thermal  expansion  of  the 
waveguide  causes  air  gaps  to  form  between  the  sample  under  test  and  the  waveguide 
walls.  This  geometry  scatters  higher-order  modes  inside  the  waveguide,  although  all 
but  the  dominant  TE10  mode  quickly  evanesce.  However,  each  higher  order  mode 
transmits  or  stores  a  finite  amount  of  complex  power  that  does  not  propagate  to  the 
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receiving  sensor  inside  the  network  analyzer.  Therefore,  characterization  using  the 
NRW  algorithm  with  the  standard  reflection  (Sn,  S22)  and  transmission  (S21,  S12) 
scattering  (S)  parameters  yields  inaccurate  results  [17,21,29].  In  addition,  the  re¬ 
flection  measurements  are  extremely  sensitive  to  the  axial  placement  of  the  sample 
in  the  waveguide.  Any  displacement  from  the  calibration  plane  along  the  waveguide 
axis  incurs  a  two-way  phase  delay  or  advance  which,  without  compensation,  yields 
inaccurate  results. 

This  thesis  presents  a  mode-matching  analysis  using  TMy  modes  to  model  a 
PFW  system  containing  a  single  air  gap  between  the  top  of  the  sample  and  the  waveg¬ 
uide.  Theoretical  S-parameters,  namely  S’ff  and  S^’,  are  calculated  and  compared 
to  the  experimental  S-parameters  obtained  from  the  network  analyzer.  A  complex 
2-D  Newton-Raphson  root  search  algorithm  solves  for  the  values  of  permittivity  and 
permeability  that  minimize  the  difference,  which  are  assumed  to  be  the  actual  values 
of  the  material  parameters.  It  is  shown  numerically  to  be  sufficient  to  enforce  con¬ 
tinuity  of  two  of  the  three  transverse  field  components  at  the  boundary  between  the 
empty  and  PFW  regions. 

To  counteract  the  error  introduced  by  shifting  the  sample  in  the  waveguide,  the 
respective  forward  and  reverse  theoretical  S-parameters  are  multiplied  together  and 
compared  to  the  experimental  forward  and  reverse  products.  A  complex  2-D  Newton- 
Raphson  root  search  algorithm  iteratively  solves  for  permittivity  and  permeability 
based  on  these  new  functions.  This  reference  plane  independent  (RPI)  formulation  is 
based  on  the  mode-matching  technique  described  above. 

The  combination  of  these  components  (mode  matching,  RPI  formulation,  and 
use  for  magnetic  materials)  has  not,  to  the  best  knowledge  of  the  author,  been  pre¬ 
sented  before. 
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1.2  Scope 

The  standard  industry  method  for  materials  characterization  uses  the  microwave 
stripline.  The  availability  of  appropriately  sized  samples,  as  well  as  ease  of  use  and 
performing  calculations,  contribute  greatly  to  this  fact.  Other  popular  test  mecha¬ 
nisms  include  the  focused  beam,  cavity  resonator,  coaxial  waveguide,  and  single-probe 
waveguide  [21].  All  of  these  methods  can  be  applied  in  high  temperature  situations, 
and  perhaps  the  focused  beam  system  gives  the  best  performance  (no  higher  order 
modes  are  excited).  Since  this  thesis  presents  a  solution  specifically  for  the  rectangu¬ 
lar  waveguide,  other  test  setups  are  not  discussed.  This  exclusion  is  appropriate,  since 
many  facilities,  including  the  sponsor  of  this  research,  do  not  possess  every  apparatus 
on  this  list. 

In  general,  materials  characterization  will  need  to  be  performed  across  a  wide 
range  of  frequencies.  However,  only  S-band  (2.6  -  3.95  GHz)  and  X-band  (8.2  - 
12.4  GHz)  measurements  are  taken  throughout  this  research.  The  mode-matching 
technique  can  easily  be  applied  to  rectangular  waveguides  at  other  frequency  bands. 
Also,  only  air  gaps  in  the  y  direction  (short  dimension)  of  the  waveguide  are  addressed; 
this  is  the  more  dramatic  problem,  since  the  transverse  TEW  E-field  is  non-zero 
at  these  gaps.  For  gaps  in  the  x  direction  (long  dimension),  the  tangential  (and 
transverse)  E-field  is  identically  zero  along  the  boundary,  so  small  gaps  have  less 
effect  on  the  measurement,  at  least  for  dielectrics  [21]. 

1 . 3  Thesis  Organization 

Chapter  2  gives  an  overview  of  previous  research  done  in  the  area  of  partially 
filled  waveguides,  and  highlights  the  contributions  of  this  work.  Chapter  3  presents 
an  review  of  electromagnetic  field  theory  necessary  for  the  research,  including  the  use 
of  vector  potentials.  The  derivation  of  the  mode-matching  solution  to  the  problem, 
as  well  as  the  inclusion  of  the  reference  plane  independent  analysis,  is  discussed  in 
Chapter  4.  Experimental  results  and  error  analysis  are  provided  in  Chapter  5,  and 
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conclusions  and  recommendations  are  summarized  in  Chapter  6. 
contain  the  NRW  algorithm  and  integral  proofs  of  Chapter  4. 


The  appendices 
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II.  Previous  Efforts 


Microwave  measurement  problems  considering  partially  filled  waveguides  have  been 
thoroughly  treated  over  the  years.  Several  texts,  such  as  those  by  Collin  [10],  Har¬ 
rington  [16]  and  Marcuvitz  [22],  address  the  issue  to  varying  degrees  of  complexity. 
This  chapter  will  review  some  of  these  methods. 

2.1  Modal  Methods 

Wexler  [28]  has  presented  a  general  analysis  of  scattered  modes  in  discontinuous 
waveguides,  without  regard  to  the  type  of  fields  present.  Boundary  conditions  and 
continuity  of  transverse  fields  are  satisfied  by  an  infinite  series  of  modes  on  each 
side  of  the  obstacle  or  junction.  As  in  this  thesis,  the  objective  is  to  determine 
the  distribution  of  complex  power  among  the  the  scattered  modes.  Wexler  uses  the 
orthogonality  relation 


on  the  non-degenerate  modes.  Although  [28]  mentions  only  PEC  obstacles,  the  modal 
formulation  may  be  used  on  dielectric  and/or  magnetic  material  obstacles  as  well. 

The  problem  of  finding  cutoff  frequencies  in  a  waveguide  partially  filled  with  an 


exponentially  varying  dielectric  has  been  addressed  by  Gonzalez  [15];  his  solution  does 


not,  however,  explicitly  solve  for  the  propagation  constants  of  higher  order  modes. 

The  problem  discussed  by  Jarern  et  al.  [19]  is  similar  to  that  treated  in  this  re¬ 
search,  although  with  notable  differences.  Using  TMy  modes,  Jarem  uses  a  method  of 
moments  analysis  to  calculate  the  theoretical  reflection  and  transmission  coefficients, 
and  uses  least-squares  curve  fitting  to  match  them  to  the  measured  S-parameter  data. 
The  axial  propagation  wave  number  7 bn  is  the  solution  of  the  eigenvalue  equation 


y/f  tanh  (7-1  nd)  =  ^  tanh  (jz2n  (d  -  b)) 

7 lin  +  7 In  +  U26biLlbi  =0,  2  =  1,2 


(2.1) 
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where  d  is  the  sample  thickness,  b  is  the  height  of  the  waveguide,  and  ^/zm  is  the 
^/-directed  wave  number  in  the  partially  filled  region.  The  index  i  —  1,2  refers  the 
subregions  of  material  (ei,/xi)  and  (€2,(12)  respectively.  To  satisfy  boundary  condi¬ 
tions  between  the  FFW  and  PFW  region,  Jarern  tries  three  methods:  matching  a 
single  E  and  H  component  at  the  FFW/PFW  interface;  using  Galerkin’s  method  to 
match  the  inner  products  of  the  fields;  and  matching  all  transverse  held  components 
individually.  Jarem  chooses  to  continue  with  the  first  method,  matching  a  single 
transverse  component  of  E  and  H.  This  choice  is  justified  by  both  numerical  cross¬ 
checks  and  a  variational  admittance  expression  which  depends  only  on  two  transverse 
components.  A  more  rigorous  validation  of  this  choice  of  boundary  conditions  is  given 
in  this  research. 

Similarly,  Catala-Civera  et  al.  [8]  present  a  method  for  extracting  complex  per¬ 
mittivity  of  a  dielectric  using  PFW  theory.  An  uncertainty  study  of  the  procedure  is 
also  given.  The  material  sample  is  discontinuous  in  x  (i.e.  long  transverse  dimension). 
An  iterative  material  perturbation  technique  is  used  to  search  for  the  correct  axial 
wave  number  7  of  the  PFW,  using 

.  /'/'  [(Ae  -  jAa)  E  •  E0  —  A/xH  •  H0]  dS 
7  -  7o  -  -jw  j j  (eo  x  H  —  E  x  H0)  •  £  dS 

where  Ae  and  A/x  are  the  material  perturbations  and  the  subscript  “0”  refers  to  the 
unperturbed  values  [6,16].  Depending  on  the  height  of  the  material  sample,  the  initial 
unperturbed  material  parameters  those  of  either  the  empty  or  full  waveguide.  The 
uncertainty  analysis  revealed  greater  error  for  low-loss  samples,  as  well  as  samples  of 
short  axial  length.  Additionally,  the  importance  of  precise  sample  alignment  with  the 
calibration  plane  is  stressed. 

Bogle  [7]  has  also  developed  a  similar  solution  to  the  partially  filled  waveguide 
characterization  problem,  although  he  uses  a  slightly  different  formulation.  The  PFW 
used  in  his  analysis  has  left-right  gaps,  with  the  sample  in  the  center.  Havrilla  has  used 
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a  perturbational  method  to  compensate  for  small  gaps  in  a  PFW,  but  has  included 
only  the  dominant  TEio  mode  in  the  solution,  excluding  higher-order  modes  [17]. 


2.2  Variational  Methods 

While  modal  methods  can  provide  an  exact  solution  for  wave  propagation  in 
an  partially  filled  waveguide,  variational  calculus  may  also  be  used  to  obtain  an  ap¬ 
proximate  solution.  Berk  [6]  outlined  a  general  variational  procedure  for  obtaining 
the  complex  propagation  coefficients  in  a  waveguide  partially  filled  with  a  dielectric 
slab.  Collin  and  Vaillancourt  [11]  have  successfully  used  the  Rayleigh-Ritz  method 
to  obtain  approximate  eigenfunctions  and  eigenvalues  in  a  waveguide  partially  filled 
with  dielectric  in  the  y- axis.  The  piecewise  function  k  (y)  corresponds  to  the  electric 
permittivity  e  as  a  function  of  position.  The  magnetic  vector  potential  is  equal  to 

A  =  fjAy  =  y  sin  ifjEn  (y)  e±lz 


where  \ pEn  (y)  is  one  of  infinitely  many  solutions  of  the  Sturm-Liouville  equation 


d^ipE  1  dn  dips 
dy 2  k  (y)  dy  dy 


+  (  Kkl  -  ^  +  R  )  IpE  =  0 


having  corresponding  eigenvalues  y2.  It  can  be  shown  that  the  solutions  ipEn  form  an 
orthogonal  set  with  respect  to  the  weighting  function  k-1.  Additionally,  the  equation 


72  L  ^/Edy-J0  I  V  dv 


d^E 


Kk0 


dy 

k{v) 


=  0,  K{y)±  0  (2.2) 


is  a  variational  expression  for  the  true  propagation  constant  y2  [2].  Collin  proposes 
the  set  of  eigenfunctions  in  the  empty  waveguide, 


,  ( nny\ 

fe>=VTcosrrJ’ 


n  —  0,1,  2... 


(2.3) 


where  e  is  Neumann  factor,  to  use  in  the  extremisation  of  (2.2).  This  method  is  used 
to  match  the  tangential  components  of  the  fields  at  the  junction  between  the  empty 
and  partially  filled  waveguide  [11].  While  this  technique  is  by  its  nature  approximate, 
it  has  the  advantage  of  avoiding  the  solution  of  a  transcendental  equation. 

The  use  of  the  Rayleigh-Ritz  method  was  improved  by  Vander  Vorst  and  Go- 
vaerts,  who  computerized  the  algorithm  for  use  in  a  variation- iteration  method  [26]. 
The  exact  solution  for  wave  propagation  in  a  waveguide  containing  .E-plane  slabs  of 
dielectric  has  been  given  by  Gardiol  [14], 

2.3  Other  Methods 

Fchlen  [13] ,  also  using  a  modal  held  expansion,  developed  a  rigorous  PFW  anal¬ 
ysis  for  samples  in  a  coaxial  test  fixture,  with  future  application  to  high-temperature 
measurements  in  mind. 

Seeking  to  improve  on  the  transmission/reflection  method  of  materials  charac¬ 
terization,  Baker-Jarvis  et  al.  developed  an  method  to  correct  errors  of  the  NRW 
algorithm  (found  in  Appendix  A)  in  low-loss  samples  with  thicknesses  approximately 
integer  multiples  of  half- wavelengths  [3] .  Relevant  to  this  thesis  is  the  presentation  of 
a  family  of  equations  that  are  independent  with  respect  to  both  the  reference  plane 
and  the  sample  thickness  itself.  These  equations  isolate  the  s-parameters  in  terms 
of  other  known  quantities,  to  be  used  in  a  minimization  equation  as  part  of  a  root 
search  for  the  correct  permittivity.  The  paper  does  not  test  magnetic  materials,  and 
is  specific  to  measurements  using  fully-filled  waveguides. 

Wilson  [29]  and  Champlin  [9]  have  analyzed  the  effect  of  an  air  gap  on  calcula¬ 
tion  of  complex  permittivity  from  transmission  and  reflection  measurements.  Wilson, 
using  Wexler’s  formulation  with  a  material  slab  discontinuity,  successfully  derived 
held  expansions  in  the  empty  and  PFW  regions.  It  is  suggested  that  a  conducting 
paste  be  applied  in  the  gap  as  a  correction,  and  this  is  done  with  great  success.  How¬ 
ever,  such  a  solution  is  not  applicable  in  a  high-temperature  PFW  situation,  since 


9 


the  paste  cannot  withstand  the  high  temperature  [21].  Champlin,  recognizing  the 
difficulty  of  obtaining  uniformity  in  the  sample,  gives  an  analytic  correction  for  small 
gaps  to  be  used  in  characterization  of  non-magnetic  materials.  He  does  not  use  modal 
analysis. 

2.4  Summary 

Notable  contributions  concerning  microwave  measurements  of  a  PFW  system 
were  reviewed  in  this  chapter.  Several  works  have  been  published  investigating  differ¬ 
ent  PFW  scenarios,  including  those  containing  variable  dielectric  material.  Solution 
techniques  are  commonly  modal  or  variational,  with  computational  results  often  given 
as  verification. 
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III.  Electromagnetics  Fundamentals 

An  understanding  of  basic  electromagnetics  and  guided  wave  theory  is  necessary  to 
develop  the  solution  for  the  problem  at  hand.  Presented  in  this  chapter  is  an  ex¬ 
planation  of  vector  potentials,  the  proper  use  as  they  pertain  to  field  construction, 
and  a  review  of  guided  wave  theory  as  it  pertains  to  this  thesis.  The  analysis  of  this 
chapter  borrows  extensively  from  [4],  [16]  and  [21].  A  eJa;t  time  dependence  is  assumed 
and  suppressed  throughout.  Readers  familiar  with  these  texts  may  feel  comfortable 
moving  to  the  next  chapter. 

3.1  Maxwell’s  Equations 

Electric  and  magnetic  field  behavior  in  simple  media,  defined  as  linear,  isotropic, 
homogeneous,  and  dispersive,  can  be  described  by  the  coupled  vector  form  of  Maxwell’s 
equations: 


V  x  E  =  — M  —  juB 

(3.1a) 

V  x  H  =  J  +  ju D 

(3.1b) 

< 

0 

S5 

(3.1c) 

V  ■  B  =  pm 

(3- Id) 

The  following  auxiliary  relations,  also  for  propagation 

in  simple  media,  apply,: 

D  =  eE  =  (e'  -  je")  E 

(3.2a) 

B  =  /iH  =  (//  —  jp!')  H 

(3.2b) 

J  =  crE 

(3.2c) 

where  u>  is  the  radian  frequency,  E  and  H  are  the  electric  and  magnetic  fields  respec¬ 
tively,  B  is  the  magnetic  flux  density,  D  is  the  electric  flux  density,  J  is  the  electric 

current  density,  M  is  the  magnetic  current  density,  pe  is  the  electric  charge  density,  pm 

is  the  magnetic  charge  density,  e  is  the  complex  electric  permittivity,  fi  is  the  complex 


11 


magnetic  permeability,  and  a  is  the  conductivity  of  the  material.  The  real  and  imag¬ 
inary  components  of  permittivity  and  permeability  are  denoted  with  single  (/)  and 
double  (//)  prime  notation,  respectively.  Relative  electric  permittivity  er  (or  magnetic 
permeability  pr)  is  defined  as  the  ratio  of  the  material  permittivity  (permeability)  to 
the  permittivity  (permeability)  of  free  space, 


e 


eo 

R 

l1  o 


Maxwell’s  equations  describe  the  coupling  of  electric  and  magnetic  fields  as  en¬ 
ergy  propagates  in  space  or  through  a  material.  The  constitutive  parameters  e,  //, 
and  a  determine  the  field  response  in  the  material  to  the  application  of  an  electro¬ 
magnetic  field  [21].  The  real  and  imaginary  parts  of  e  and  /i  represent  energy  stored 
in  the  material  and  loss  mechanisms  of  the  material,  respectively. 

3.2  Vector  Potentials 

It  is  useful  to  define  auxiliary  functions  to  aid  in  the  solutions  of  problems 
involving  Maxwell’s  equations,  such  as  the  partially  filled  waveguide  problem  of  this 
thesis. 

3.2.1  Magnetic  Vector  Potential.  The  fields  generated  by  an  electric  current 
in  a  region  free  of  magnetic  sources  (i.e.  M  —  0)  must  satisfy  Gauss’s  Law,  V  ■  B  =  0. 
Since  in  a  source  free  region  B  is  always  solenoidal,  this  implies  it  has  a  magnetic 
vector  potential  A,  such  that 


B  4  =  V  x  A 


(3.3) 


or,  by  (3.2b), 


(3.4) 
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Substituting  (3.4)  into  (3.1a)  yields 


V  x  E_4  =  =  —  jwV  x  A 

which  can  also  be  written  as 

V  x  [Ej4  +  ju A]  =  0 

By  identity,  this  implies  that  E^  +  jujA  has  a  scalar  potential,  that  is 

Ea  +  juA  =  -V0e 


(3.5) 


(3.6) 


(3.7) 


or, 

E  a  =  —  V(/>e  —  jtuA  (3.8) 

where  the  electric  scalar  potential  (f>e  is  a  function  of  position. 

By  taking  the  curl  of  both  sides  of  (3.4)  and  using  the  vector  identity  V  x  V  x 
A  =  V  (V  •  A)  -  V2A,  it  can  be  reduced  to 

/iVx  H  =  V  (V  •  A)  -  V2A  (3.9) 

Using  (3.1b)  leads  to 

fiJ  +  jujefiEiA  =  V  (V  •  A)  —  V2A  (3.10) 


Substituting  (3.8)  into  (3.10)  obtains 

V2A  +  k2 A  =  —fi J  +  V  (V  ■  A)  +  V  (jujefj,(l>e)  (3.11) 

=  —  fiJ  +  V  ( V  ■  A  +  jujt lief) e)  (3-12) 
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where  k2  =  a;2e^.  Having  previously  defined  the  curl  of  A  in  (3.3),  the  divergence  of 
A,  which  is  independent  of  the  curl,  may  be  defined  as 

V  •  A  =  —jujefi(j)e  (3.13) 

which  is  known  as  the  Lorentz  gauge  condition.  Substituting  into  (3.11)  results  in 

V2A  +  k2  A  =  -/jJ  (3.14) 


This  is  known  as  the  wave  equation  for  A.  The  electric  vector  field  due  to  the  magnetic 
vector  potential  A,  using  (3.13)  in  (3.8),  can  be  written  as 

=  —  V0e  —  joj A  =  —ju A  H - V  (V  •  A)  (3.15) 

jtvefi 

Thus  knowledge  of  A  enables  the  finding  of  and  from  (3.15)  and  (3.4),  re¬ 
spectively.  The  components  of  the  fields  due  to  an  electric  current  density  J  having 
been  found,  it  remains  necessary  to  calculate  the  fields  due  to  a  magnetic  current  M. 


3. 2. 2  Electric  Vector  Potential.  The  fields  generated  by  an  equivalent  mag¬ 
netic  current  in  a  region  free  of  electric  sources  (i.e.  J  =  0)  must  satisfy  Gauss’s  Law, 
V  •  D  =  0.  Observing  that  D  is  also  solenoidal,  this  implies  it  has  an  electric  vector 
potential  F  such  that 

DF  =  V  x  F  (3.16) 


or,  by  (3.2a), 


E/7  = 


— -V  x  F 

e 


Substituting  (3.17)  into  (3.1b)  yields 


(3.17) 


V  x  Hf  =  jujcEp  =  ju\7  x  F 


(3.18) 
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which  can  also  be  written 

V  x  [Hf  +  jujF]  =  0  (3.19) 

By  identity,  this  implies  that  +  juF  has  a  scalar  potential,  that  is 

HF  +  juF  =  -V(j)m  (3.20) 

or 

Hp  =  —  Vcftm  —  jwF  (3-21) 

where  the  magnetic  scalar  potential  is  a  function  of  position. 

By  taking  the  curl  of  both  sides  of  (3.17)  and  using  the  vector  identity  V  x  V  x 
F  =  V  (V  •  F)  —  V2F,  it  can  be  reduced  to 

-eV  x  Ef  =  [V  (V  ■  F)  -  V2F]  (3.22) 

Using  (3.1a)  leads  to 

eM  +  jue/iFlp  —  V  (V  •  F)  —  V2F  (3.23) 

Substituting  (3.21)  into  (3.23)  obtains 

V2F  +  k2 F  =  — eM  +  V  (V  ■  F  +  juefMfrm)  (3.24) 

Having  previously  defined  the  curl  of  F  in  (3.16),  the  divergence  of  F,  which  is 
independent  of  the  curl,  may  be  defined  as 

V  •  F  =  -jijje/Mpm  (3.25) 

which  is  known  as  the  Lorentz  gauge  condition.  Substituting  this  into  (3.24)  results 
in 

V2F  +  k2  F  =  -eM  (3.26) 
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This  is  the  wave  equation  for  F.  The  magnetic  vector  field  due  to  the  electric  vector 
potential  F,  using  (3.25)  in  (3.21)  can  be  written  as 

Hp  =  — - — V  (V  •  F)  —  juF  (3.27) 

juiefi 

Thus  knowledge  of  F  enables  the  determination  of  Ep  and  Hp  from  (3.17)  and  (3.27), 
respectively. 

3.2.3  Summary.  It  has  been  shown  that  the  electromagnetic  field  compo¬ 
nents  sustained  by  an  electric  current  density  J  or  magnetic  current  density  M  can  be 
calculated  through  use  of  the  magnetic  vector  potential  A  or  electric  vector  potential 
F.  When  both  sources  are  present,  the  principle  of  superposition  may  be  applied  to 
determine  the  total  fields,  namely 


E  total  —  E^  +  Ep  (3.28) 

and 

H  total  =  H  4  +  Hp  (3.29) 

Using  vector  potentials  to  construct  solutions  to  Maxwell’s  equations  is  the 
subject  of  the  next  section. 

3.3  Field  Construction 

The  simplest  solutions  to  Maxwell’s  equations  are  Transverse  ElectroMagnetic 
(TEM)  field  configurations,  or  modes,  where  both  the  electric  and  magnetic  field 
components  are  transverse  to  the  direction  of  propagation.  Modes  that  are  Transverse 
Electric  (TE)  or  Transverse  Magnetic  (TM)  may  also  be  constructed.  The  direction 
of  propagation  is  indicated  by  superscript,  i.e.,  TEM2  is  transverse  to  the  z  direction. 
In  this  chapter,  only  solutions  in  a  rectangular  coordinate  system  are  considered. 
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3.3.1  TEM  Modes.  In  a  source  free  region,  both  J  and  M  do  not  exist. 
Equations  (3.14)  and  (3.26)  become  the  homogeneous  differential  equations 

V2A  +  k2  A  =  0  (3.30) 

and 

V2F  +  k2  F  =  0  (3.31) 

Since  (3.30)  and  (3.31)  are  of  the  same  form,  the  solutions  will  also  be  of  the  same 
form.  Therefore,  the  following  development  is  only  be  presented  for  A.  Letting 


A  (x,  y,  z)  =  xAx  (x,  y,  z )  +  yAy  (x,  y,  z)  +  zAz  (x,  y,  z)  (3.32) 

the  component  form  of  (3.30)  is 


V2AX  +  k2Ax  =  0 

(3.33a) 

V2Ay  +  k2Ay  =  0 

(3.33b) 

V2A,  +  k2Az  =  0 

(3.33c) 

The  total  electric  held  is  given  by  the  sum  of  (3.15)  and  (3.17), 


Etotal  —  E^4  +  Ep  —  —ju  A 


— — V  (V  ■  A)  -  -V  x  F 

juey  e 


(3.34) 


Applying  (3.32)  to  (3.34),  the  total  E  held  can  be  written  in  component  form 


E  = 


x 


y 


-juAx  +  -A_  +  Uf) 

jusiy  -r  juiey  ^  dxdy 


,  ,  PaA 

-r  a..2  -T  dydz  j 


dF, 

dy 


dFy 

dz 


dy 2 


1  (dF*  _  9F\ 


V  dz 


dx 


'  -jujAz  +  +  + 

J  A  juey  \  oxoz  oyoz  ozz  I  e  \  ox  oy  J 


(3.35) 
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By  the  Duality  Principle  the  total  magnetic  field  likewise  is  given  by 


Htotal  =  Ha  +  HF  =  +  —V  (V  •  F)  -  jujF 

which  which  can  also  be  expanded  in  components  as 


H  =  x 
y 

z 


-jujFx  + 


1  ( 

jcoe^L  l  dx 2 


-jUjFy  + 

~h  jujey  V  dxdz 


(  + 


PFy 
'  dxdy 

d2Fl 
dy 2 

d2F \ 
dydz 


,  d^FA  |  1  (dA, 

'  dxdz  J  ‘  fi  y  dy 


y~  + 
y  + 


d2Fz 

dydz 

d2Fz 

dz2 


+ 


I  (9AX 

y,  V  dz 
1  ( 9Ay 


dx 


dAz\ 
dx  J 


(3.36) 


(3.37) 


To  generate  a  TEMZ  mode,  both  A  and  F  must  be  used  as  generating  functions. 
Requiring  A  =  zAz  and  F  =  zFz  to  be  non- zero  while  Ax  =  Ay  =  Fx  =  Fy  =  0  will, 
by  (3.34)  and  (3.36),  produce  the  required  TEM2  field  components. 


3.3.2  TM,TE  modes.  While  the  use  of  TEM  modes  is  sufficient  for  many 
applications  (such  as  wave  propagation  in  free  space),  other  boundary  conditions 
require  the  use  of  TM  and  TE  modes.  This  nomenclature  indicates  that  either  the 
magnetic  or  electric  field  components  lie  in  a  plane  transverse  to  a  given  direction. 
For  example,  a  TEy  field  configuration  implies  that  Ey  =  0;  the  remaining  electric 
field  components,  and  all  the  magnetic  field  components,  may  or  may  not  exist. 

To  generate  modes  TM  (or  TE)  to  a  given  direction,  it  is  sufficient  to  allow  the 
magnetic  vector  potential  A  (or  electric  vector  potential  F)  to  have  a  single  non-zero 
component  in  the  direction  in  which  the  fields  are  desired  to  be  transverse  [4].  For 
example,  the  generating  potential  for  TM2  field  components  is 


A  =  (x,  y,  z) 

F  =  0 


(3.38) 
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where  the  vector  potential  A  must  satisfy  the  wave  equation  (3.14)  with  J  =  0.  Using 
the  separation  of  variables  technique,  solution  is  assumed  to  be  in  the  form 

Az  (x,  y,z)  —  f  (x)  g  (y)  h  (2)  (3.39) 

where  the  functions  f,g  and  h  must  be  chosen  to  satisfy  the  wave  equation  (3.33c) 
and  boundary  conditions  of  the  problem.  For  a  rectangular  waveguide  oriented  on 
the  ^-axis  it  is  easiest  to  apply  the  boundary  conditions  if  the  solution  is  written  in 
the  form 

Az  =  [Ci  cos  (kxx)  +  Di  sin  (kxx)]  [C2  cos  ( kyy )  +  D2  sin  (kyy)]  (yB+e~lzZ  +  B~elzZ ) 

(3.40) 

where  the  complex  exponentials  represent  traveling  waves  and  the  sine  and  cosine 
functions  represent  standing  waves.  The  separation  of  variables  also  obtains  the  con¬ 
straint  equation 

kl  +  k2y  -  7,2  =  k2  =  A/z  (3.41) 

In  (3.40),  the  coefficients  C ,  D ,  and  B  are  amplitude  constants  determined  upon  satis¬ 
faction  of  the  boundary  conditions.  The  wave  numbers  kx  and  ky  are  spatial  constants 
that  describe  field  variation  in  x  and  y  respectively,  while  the  wavenumber  k  is  de¬ 
pendent  on  both  frequency  and  the  fundamental  properties  of  the  medium  itself.  The 
wavenumber  7^  describes  held  behavior  in  the  propagation  axis,  and  consists  of  a 
real  and  imaginary  part,  such  that  7  =  a  +  j/3.  The  real  part  a  is  the  attenuation 
constant,  and  the  imaginary  part  (5  is  the  phase  constant. 

The  transverse  wave  numbers  may  also  be  complex,  and  are  often  referred  to 
as  eigenvalues,  since  they  are  characteristic  solutions  to  the  wave  equation.  The 
notation  e~lzZ  indicates  a  wave  traveling  in  the  forward,  +2,  direction,  while  elzZ 
indicates  reverse  travel  along  —2.  If  is  purely  imaginary,  the  wave  is  unattenuated. 
If  7 z  is  purely  real,  the  wave  is  evanescent.  For  complex  7-  and  Re  (72)  >  0,  the  wave 
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Figure  3.1:  Coordinate  system  used  thoughout  this  thesis.  The  waveguide  axis  is 

z,  y  is  top  to  bottom,  and  x  is  left  to  right.  The  waveguide  is  filled  with  free  space, 
with  parameters  (eo,ho) 

travels  within  some  envelope  of  attenuation.  Once  Az  is  found,  the  E  and  H  field 
components  can  be  determined  according  to  (3.35)  and  (3.37). 

3.4  Guided  Waves 

This  section  presents  the  method  for  solving  Maxwell’s  equations  in  a  rectan¬ 
gular  hollow  waveguide  with  (assumed)  perfect  electric  conducting  (PEC)  walls,  and 
cross  section  uniform  to  the  direction  of  propagation.  The  analysis  throughout  this 
section  is  similar  to  Chapter  8  of  [16].  The  geometry  of  the  waveguide  is  given  in 
Figure  3.1.  Let  the  axis  of  the  waveguide  be  the  z  axis,  with  (x,y)  dimensions  from 
the  origin  to  (a,  5);  usually,  a  ~  2b,  where  b  is  the  height  of  the  waveguide.  The 
interior  of  the  waveguide  will  be  free  space  with  electrical  parameters  60  and  [Xq.  A 
TEM  mode  cannot  propagate  in  a  hollow  waveguide  due  to  the  absence  of  both  the 
axial  currents  and  magnetic  flux  necessary  to  generate  transverse  field  components. 
Therefore,  solutions  TM  and  TE  will  be  investigated. 

Transverse  magnetic  to  2  modes,  as  previously  mentioned,  are  generated  from 
A  =  zAz  and  F  =  0.  Standing  waves  will  occupy  the  transverse  dimensions,  while  a 
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traveling  wave  will  exist  along  z.  Therefore,  use  of  the  magnetic  vector  potential  in 
(3.40)  is  appropriate.  Boundary  conditions  require  that  the  tangential  component  of 
the  electric  field  (Ez)  vanish  at  x  =  0,  x  =  a,  y  =  0  and  y  =  b.  Therefore,  /  and  g  of 
(3.39)  are 

/  (x)  =  D1  sin  (kxx)  kx  =  ™  m  =  1,2,3... 
g  (y)  =  D2  Sin  (kyy)  ky  =  n  =  1,2,3... 

Each  pairing  of  the  integers  m  and  n  represent  a  possible  mode.  If  only  forward- 
traveling  waves  are  considered,  the  TM2  mode  functions  are 


Az  =  Knn  sin 


/ mux 
V  a 


sin 


”^1)  e-7.* 


(3.42) 


where  the  constants  have  been  combined.  The  constraint  equation  (3.41)  then  be¬ 
comes 

(v)2+(t  )2-X-=*0=^.  (3.43) 

The  field  components  are  obtained  from  (3.34)  and  (3.36).  A  solution  using  TE  modes 
can  be  determined  in  a  similar  manner. 


3. 4-1  Impedance.  The  use  of  the  term  impedance  to  describe  the  complex 
ratio  of  voltage  and  current  was  applied,  in  turn,  to  circuit  theory,  transmission  lines, 
and  electromagnetic  fields  [24],  Several  representations  of  impedance  exist,  depending 
on  the  field  type,  direction  of  travel,  and  the  medium  of  propagation.  The  intrinsic 
impedance  rj  of  a  material,  dependent  only  the  parameters  e  and  p,  is 


V  = 


(3.44) 


The  wave  impedance, 


-I 


(3.45) 


is  the  ratio  of  transverse  components  of  E  to  transverse  components  of  H.  The 
impedance  is  dependent  on  the  type  of  wave:  TEM,  TE,  TM  all  have  characteris- 
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tic  wave  impedances  particular  to  the  type  of  waveguide,  material,  and  frequency.  If 
a  plane  wave  (TEM,  TE,  or  TM)  is  normally  incident  on  the  material,  then  the  wave 
impedance  Z  is  reduces  to  the  intrinsic  impedance  77. 


3.4.2  Cutoff  Frequency.  Forward  wave  propagation  (without  attenuation) 
in  the  waveguide  occurs  when  qz  =  j/3  is  purely  imaginary.  When  qz  is  purely  real, 
the  wave  is  evanescent  and  does  not  propagate.  The  transition  between  these  two 
states  occurs  at  the  cutoff  frequency,  when 


k2  ~  (kl  +  K)  =k2  -k2c  =  0 


(3.46) 


Since  k  =  =  27r/v/e/i,  substituting  in  the  wavenumbers  from  (3.42)  and 

solving  (3.46)  for  frequency  /  results  in 


(3.47) 


which  is  an  expression  for  the  cutoff  frequency  of  a  TMmn  mode.  Modes  (TM,  TE,  or 
both)  that  have  the  same  cutoff  frequency  are  called  degenerate.  The  mode  with  the 
lowest  cutoff  frequency  in  a  particular  guide  is  the  dominant  mode.  In  a  rectangular 
waveguide  where  a  >  b,  the  TEf0  mode  is  dominant. 

3-4-3  Mode  Orthogonality.  The  solution  to  the  PFW  problem  requires  the 
use  of  an  arbitrary  number  of  held  modes.  The  total  held  in  the  waveguide  is  the 
superposition  of  all  the  modes,  and  each  mode  vector  e*  is  orthogonal  to  all  other 
mode  vectors.  If  two  modes,  either  TE  or  TM,  are  considered,  multiplying  them 
together  forms  the  product 


ei  •  ej  =  Vt*i  ■ 


(3.48) 
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where  $  is  a  mode  function  and  Vt  is  the  transverse  gradient.  By  Green’s  first 
identity, 


'f'i'&jds 


(3.49) 


Also,  by  Green’s  second  identity, 


[(feci) 


2 


'f'i'&jds 


0 


(3.50) 


If  kci  7^  kCj  the  integral  must  vanish,  as  must  the  right  hand  side  of  (3.49),  that  is 


■  ej  ds 


o,  « /  :i 


(3.51) 


It  is  equivalent  to  state  that  if  the  inner  product  of  two  functions  (or  vectors)  is  zero, 
then  they  are  orthogonal.  The  proof  (from  Chapter  8  of  [16])  may  be  extended  to  the 
magnetic  held  vectors  h  as  well.  The  transverse  electromagnetic  held  at  any  point 
can  then  be  expressed  as  the  sum  of  the  mode  vectors: 


E  t  = 

i 

H,  =  £h, 

i 

where  the  mode  vectors  e  and  h  can  be  either  TE  or  TM. 


3-4-4  Power  Transmission.  When  many  modes  in  a  waveguide  with  PEC 
walls  exist  simultaneously,  each  will  propagate  energy  independently  [16].  For  the 
current  test  setup  of  a  PFW  this  is  especially  relevant,  since  it  is  expected  that  an 
infinite  number  of  TMy  modes  will  exist  ou  both  sides  of  the  discontinuity.  In  the 
fully  filled  waveguide,  all  the  energy  is  transmitted  by  the  dominant  TE10  mode;  this 
allows  the  use  of  closed  form  algorithms,  such  as  NRW,  which  consider  single-mode 
propagation.  When  multiple  modes  exist,  it  is  necessary  to  take  into  account  the 
complex  power  transmitted  by  both  the  propagating  and  evanescent  modes.  For  N 
modes,  the  total  ^-directed  complex  power  in  a  section  of  waveguide  is  the  sum  of  the 
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power  contained  in  each  mode,  such  that 


=  E  v‘ri  <3-52) 

i 

where  Vi  is  the  mode  voltage  and  is  the  mode  current  of  the  ith  mode.  It  is 
apparent  that  this  is  the  expected  result,  since  the  complex  power  of  any  system  is 
simply  Pz  —  V  ■  I* . 

3. 5  Summary 

A  cursory  review  of  electromagnetics  was  presented  in  this  chapter,  including 
the  use  of  vector  potentials  and  the  modal  expansion  of  fields  in  a  waveguide.  These 
concepts  form  the  basis  of  the  modal  analysis  found  in  Chapter  IV. 
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IV.  Mode-Matching  Analysis 

The  characterization  of  electromagnetic  materials  by  transmission  and  reflection  meth¬ 
ods  in  a  rectangular  waveguide  assumes  the  sample  to  be  homogeneous  and  precisely 
machined,  i.e.  it  completely  fills  the  inner  dimensions  of  the  waveguide  and  is  perpen¬ 
dicularly  planar  to  the  waveguide  walls,  as  seen  in  Figure  4.1.  At  high  temperatures, 
thermal  expansion  distorts  the  geometry  of  both  the  sample  and  the  waveguide  it¬ 
self  [21],  The  metal  waveguide  tends  to  expand  at  a  greater  rate  than  the  sample 
under  test,  resulting  in  air  gaps  between  the  sample  and  the  waveguide  walls.  Figure 
4.2  (a)  is  an  example  of  a  top-bottom  gap  in  the  short  dimension  of  the  waveguide. 
The  boundary  conditions  at  the  gap  require  the  excitation  of  higher  order  modes 
which  are  not  accounted  for  in  closed  form  solutions  such  as  the  NRW  algorithm, 
which  uses  the  dominant  (TE10)  mode  only  [17,21], 

Development  of  the  solution  for  a  single  top  gap  will  follow  that  suggested  by 
Collin  [10]  and  Harrington  [16],  using  TMy  modes  to  construct  the  electric  and  mag¬ 
netic  fields  in  the  unobstructed  and  partially  filled  waveguide  regions.  By  applying 
appropriate  boundary  conditions,  and  ensuring  that  tangential  components  of  E  and 
H  are  matched,  a  system  of  equations  that  accurately  describes  the  structure  can 
be  developed.  In  turn,  the  forward  reflection  ( S ^y)  and  transmission  (S^y)  coeffi¬ 
cients  can  be  extracted.  A  minimization  of  the  difference  between  the  theoretical  and 
experimental  S-parameters  is  calculated  according  to  the  minimization  equations 


S'^y  (u,  e,  fi)  —  S^p  (u>)  <  tol 
Sfi  (w,  e,  /i)  -  ^2iP  M  <  to 1 


(4.1) 


A  Newton-Raphson  root  search  of  the  complex  parameters  e  and  fi  is  performed  until 
an  acceptable  minimum  tolerance  is  reached.  In  addition,  it  is  necessary  to  have 
available  a  reference  plane  independent  measurement  scheme  to  compensate  for  error 
in  sample  placement  or  if  alignment  with  the  calibration  plane  cannot  be  guaranteed. 
This  method  is  presented  in  the  final  section  of  the  chapter. 


25 


Figure  4.1:  Forward  scattering  parameters  in  a  fully  filled  waveguide  measurement. 
The  sample  is  aligned  with  calibration  plane  1;  calibration  plane  2  is  the  other  end  of 
the  sample  holder. 

4-1  Newton- Raphson  Root  Search 

The  minimization  equations  used  in  this  thesis  can  be  thought  of  as  two  func¬ 
tions,  dependent  on  values  of  er  and  fir  that  force  the  differences  in  (4.1)  below  the 
required  tolerance.  The  complexity  of  the  desired  root  search  depends  on  the  parame¬ 
ters  of  the  material.  If  a  material  has  both  electric  and  magnetic  loss  properties,  then 
a  2-dimensional  root  search,  using  two  equations  and  two  unknowns,  must  be  used. 
However,  if  a  material  is  a  non-magnetic  dielectric,  only  er  is  the  unknown  (pr  =  1), 
and  a  complex  1-dimensional  root  search  of  one  equation  and  one  unknown  can  be 
applied.  This  technique  will  be  shown  first. 

One  of  the  most  common  techniques  in  numerical  analysis  used  to  determine 
roots  of  equations  is  the  Newton-Raphson  method  [5].  If  p  is  an  unknown  root,  and  a 
function  /  is  differentiable  on  the  interval  of  all  approximations  to  p,  then  f  (p)  =  0. 
Let  x  =  Xq  be  an  initial  guess  of  the  root  p.  A  Taylor  series  expansion  of  /  (x)  around 
the  approximate  root  Xo  is 

At2 

f(x)  =  f  ( x0  +  A x)pa  f  (x0)  +  A xf  (x0)  +  (0  +  ...  (4.2) 
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Cal  Plane  1 


Cal  Plane  2 


(b) 


Figure  4.2:  Two  cross  sectional  views,  (a)  transverse  and  (b)  axial,  of  a  partially 

filled  waveguide  consisting  of  a  single  top  air  gap,  with  material  parameters  (e,p). 
Only  the  dominant  mode  forward  scattering  parameters,  Su  and  S21  are  shown,  al¬ 
though  the  PFW  geometry  excites  infinitely  many  modes. 
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where  £  is  on  the  range  x0  to  x0  +  Ax.  If  the  approximation  x0  +  Ax  is  set  equal  to 
p,  then 

Ax2 

f  (p)  =  0  «  /  (x0)  +  A xf  (x0)  +  -7 jp/"  (0  +  (4.3) 

he  hrst  derivative  of  /  can  be  calculated  numerically,  using  forward,  central,  or  back¬ 
ward  differences.  If  Ax  is  small,  then  the  third  and  higher  order  terms  can  be  ne¬ 
glected.  This  leads  to  an  estimate  for  Ax, 


Ax 


f(z  0) 

/'  (x0) 


(4.4) 


The  next  approximation  Xi  to  the  root  is  obtained  by  adding  Ax  to  the  previous 


estimate,  xq, 


x\  =  xq  +  Ax  =  xq  — 


f(z  0) 

f  (zo) 


In  general,  then,  the  nth  approximation  to  the  root  p  is 


(4.5) 


xn 


%n—  1 


f{Zn- 1) 
/'  (®n-l) 


(4.6) 


Convergence  is  dependent  on  a  “good”  initial  guess  for  the  root  p.  If  the  root  is  known 
to  be  complex,  the  initial  guess  must  also  be  complex;  the  method  will  not  converge 
to  a  complex  root  if  a  purely  real  initial  guess  is  supplied  [1]. 

If  roots  to  more  than  one  equation  are  to  be  found  simultaneously,  Newton- 
Raphson’s  method  can  easily  be  expanded  to  higher  dimensions.  Consider  two  func¬ 
tions 


f(u,v)  =  0 

g{u,v)  =  0  (4.7) 


with  roots  u  and  v.  It  is  assumed  both  functions  are  differentiable  on  the  interval,  or 
surface,  of  approximations  to  u  and  v.  If  x  =  xq  and  y  =  yo  are  supplied  as  estimates 


to  the  roots,  then  the  Taylor  series  expansion  around  the  approximations  is 


d  f  df 

f  (x,  y)  =  f  {xq  +  Ax,  y0  +  Ay)  «  /  (x0,  y0)  +  Ax  +  —  Ay 

ox  ay 

9  (x,  y)  =  g  (x0  +  Ax,  y0  +  A y)  ^  g  (x0,  y0)  +  ||  Ax  +  ^ Ay  (4.8) 

where  the  higher  order  derivatives  of  the  series  have  been  neglected.  Letting  xo+ Ax  = 
u  and  y0  +  Ay  =  v,  (4.8)  reduces  to 


df  d  f 

f  (u,  v)  =  0  =  /  (x0  +  Ax,  y0  +  Ay)  =  f  (x0,  y0)  +  /Ax  +  -^-A y 

dx  dy 

da  da 

g  (u,  v)  =  0  =  g  (x0  +  Ax,  y0  +  Ay)  =  g  (x0,  y0)  +  ^Ax  +  —  Ay  (4.9) 


which  can  be  recast  in  matrix  notation  as 


fx  fy 

Ax 

f(x0,y0) 

9x  9y 

Ay 

g(xo,yo ) 

The  matrix  can  be  inverted  to  obtain  expressions  for  Ax  and  Ay, 


(4.10) 


Ax 

-1 

9y 

~fv 

f(xo,yo) 

Ay 

1 

9x 

fx 

g{x0,yo) 

(4.11) 


Proceeding  as  before,  the  next  set  of  estimates  to  the  roots  are  found  by  adding  Ax 
and  Ay  to  the  previous  estimates,  xq  and  y q. 


X\  =  xq  T  Ax 

yi=yo  +  &y 


(4.12) 


The  nth  root  approximation  is  therefore 


% n  %n—  1  T  ^•%n— 1 
Vn  Vn—l  T  ^Vn— 1 


(4.13) 
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Depending  on  the  slope  of  each  function  /  and  g,  the  initial  guesses  usually  must  be 
good  approximations  to  the  root.  If  the  third  term  of  the  Taylor  series  is  too  large, 
then  it  may  be  necessary  to  include  it  in  order  to  achieve  root  convergence  [5]. 


4-2  Partially  Filled  Waveguide  Technique 

4-2.1  TAP  Modes.  In  rectangular  waveguides  where  a  ~  2b,  the  dominant 
mode  is  TE10;  the  electric  field  component  is  oriented  in  the  y  direction,  and  both 
the  x  and  z  components  do  not  exist.  Conversely,  the  magnetic  field  has  non-zero 
components  in  the  x  and  z  directions,  while  the  y  component  does  not  exist. 

Recall  from  [16]  that  the  components  of  the  E-field  of  the  dominant  TE10  mode 
can  be  written  as 

ex  0 

ey  =  E0  sin  ( kxx )  e~lzZ 
ez  =  0 

and  the  H-field  components  as 


hx  =  jj^Eq  sin  ( kxx )  e~lzZ 

hy  =  0 

hz  =  ~^E0  cos  ( kxx )  e~lzZ 


The  wave  impedance,  by  definition,  is  the  ratio  of  transverse  e  and  h  vector  field 
components  of  the  dominant  mode  and  is  therefore  equal  to 


7  _  y 

ZjTEz  —  7T 

10  ht 


Eq  sin  (kxx)  eTJzZ  jeoy 0 

^E0  sin  (kxx)  e^*z  7 * 


(4.14) 


Since  this  is  equal  to  the  wave  impedance  of  a  TM^0  mode,  and  because  all  the 
scattered  modes  will  be  in  the  set  of  TMfn,  the  complete  mode  set  can  be  constructed 
using  TMy  modes. 
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Figure  4.3:  The  three  computational  regions  used  in  held  construction:  Region  I: 
empty;  Region  II:  PFW;  Region  III:  empty. 

4-2.2  Field  Construction.  It  is  sufficient  to  use  the  magnetic  vector  potential 
A  to  generate  the  necessary  electromagnetic  held  components  required  for  satisfying 
boundary  conditions  in  all  the  regions  of  the  waveguide.  Figure  4.3  is  a  diagram  of 
the  regions  discussed  in  this  section.  In  the  two  empty  regions  (i.e.  Regions  I  and 
III),  the  magnetic  vector  potential  takes  the  form 

Ay  =  ( Bi  sin  kxx  +  B2  cos  kxx )  ( C\  sin  kyy  +  C2  cos  kyy )  (.D+e~702  +  D“e702)  (4.15) 

where  kx,  ky ,  and  70  are  the  transverse  wave  numbers  and  complex  propagation  wave 
number,  respectively,  that  satisfy  the  constraint  equation  in  free  space, 

kl  +  kl-ll  =  kl  = 

The  use  of  e_7°2  denotes  a  forward  propagating  wave  in  the  z  direction  (i.e.  through 
the  waveguide),  and  70  is,  in  general,  a  complex  quantity.  The  respective  E  and 
H-helds  can  be  generated  by 

Ea  =  -7 — - —  [V  (V  •  A)  +  /cqA]  ,  =  ^  (4.i6) 

J^eo/Uo  yo 
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Substituting  (4.15)  into  (4.16)  and  simplifying,  the  E  components  may  be  writ¬ 
ten  as 


e*  = 


{Bi  cos  kxx  -  B2  sin  kxx)  (Ci  cos  kyy  -  C2  sin  kyy )  {D~e~loz  +  D~eloz ) 

=  ^ kj^c^]  (Bi  sin  kxx  +  B2  cos  kxx)  (Ci  sin  kyy  +  C2  cos  kyy )  ( D+e~loZ  +  D~eloZ) 

(Bi  sin  kxx  +  B2  cos  kxx)  (Ci  cos  kyy  —  C2  sin  kyy )  (D+e_7°2  —  D~eloZ) 

(4.17a) 


y  jueoW) 

 -ky'ro 


jojeofio 


and  the  H  components  as 


hx  =  ~  (Bi  sin  kxx  +  B2  cos  kxx)  (C\  sin  kyy  +  C2  cos  kyy)  (D+e~loz  —  D~e'roz ) 
hy  =  0 

hz  =  (B\  cos  kxx  —  B2  sin  kxx)  (C\  sin  kyy  +  C2  cos  kyy )  (D+e_7°2  +  D~eloZ) 

(4.17b) 

Satisfaction  of  electric  held  boundary  conditions  at  the  waveguide  walls  requires  B2 
and  C i  to  vanish,  leading  to  the  definition  of  the  wave  numbers  kx  and  ky : 


l,  __  "wr  u  n-n 

—  a  >  Fvy  b 

m  =  1,  2,  3,  ... 
n  =  0, 1,  2, ... 

where  a  is  the  long  dimension  of  the  waveguide  cross  section,  and  b  is  the  short 
dimension,  x  and  y  respectively  (refer  to  Figure  4.2  (a)).  By  the  Uniqueness  Theorem, 
it  is  sufficient  to  satisfy  only  the  electric  held  condition  [4], 

The  complex  mode  propagation  wave  number  y0  of  the  empty  region  of  the 
waveguide  can  then  be  written  as 
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By  applying  the  aforementioned  boundary  conditions,  the  Region  I  electric  and  mag¬ 
netic  field  components  of  (4.17)  can  simplified  as 

ex  =  ~ue0kJ0  cos  kxx  sin  kyy  ( D+e~loZ  +  D~eloZ ) 
ey  =  oy]  s^n  ^xx  cos  kyy  ( D+ e~loZ  +  D~eyoZ) 

_  ky 70  s[n  x s[n  k  y  ( D+e~loz  —  D~eyoz) 

z  J.  yv  \  j 

hx  =  sin  kxx  cos  kyy  ( D+e~loZ  —  D~eloZ) 

hy  =  0 

hz  —  —  cos  kxx  cos  kyy  ( D+e~loz  +  D~eyoz ) 


Using  PFW  theory  in  Region  II,  spatial  shift  factors  are  used  to  create  an 
alternative  magnetic  vector  potential  A,  which  is  then  applied  to  construct  the  fields 
in  each  subregion  (material  and  free  space),  namely 


Ayi  =  sin  kxx  ( Cx  sin  kyly  +  Dx  cos  ky±y )  (R+e  lnZ  +  Bn  elnZ ) 

(4.19) 

Ay2  =  sin  kxx  {C2  sin  ky2  (b  -  y)  +  D2  cos  ky2  ( b  -  y))  (C+e~7nZ  +  C~e7nZ ) 

where  the  subscripts  1  and  2  represent  each  subregion. 

In  Region  II,  the  electric  field  boundary  condition  in  x  is  identical  to  that  of 
Region  I,  namely,  the  tangential  fields  vanish  at  x  =  0  and  x  =  a.  Therefore,  a 
standing  sine  wave  exists  in  both  subregions  of  Region  II  and  throughout  Region  I; 
this  is  manifested  in  the  choice  of  harmonic  functions  in  (4.19).  From  (4.16)  the  fields 
in  each  region  can  be  calculated.  If  the  sample  has  height  h,  then  fields  in  the  first 
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subregion  (0  <  y  <  h)  exist  within  the  material: 


exi  =  — cos  kxx  (Ci  cos  kyly  -  Dx  sin  kyly)  (B+e  7n2  +  Bn  e7"2) 

jueiHi 
(fr  2  _  £  2  1 

eyi  =  — 7 - —  sin  kxx  (Ci  sin  kyiy  +  Di  cos  kyly)  (B+e_7ri2  +  5“e7"2) 

Ci  =  -t — — —  sin  kxx  (Ci  cos  kyiy  —  Di  sin  ky\y)  (B~ j~ e-7"2  —  S“e7"2) 

jujeiyi 

hx i  =  —  sin  /cxx  (Ci  sin  kyly  +  Efi  cos  /cyiy)  (B+e“7"2  -  B~elnZ ) 
hi 

/iyi  0 

k 

hz\  =  —  cos  kxx  (Ci  sin  kyly  +  Z?i  cos  /q^y)  (-B^e_7nZ  +  B~elnZ^  (4.20) 

hi 

Fields  in  the  second  subregion  (h  <  y  <b)  exist  in  free  space: 

ex2  =  .  k'l  ky2  cos  kxx  (C2  cos  ky2  (b  -  y)  -  D2  sin  ky2  (b  -  y))  (C+e_7n2  +  C~elnZ ) 
joieoho 

(tf  _  p  i 

ey2  =  — - - —  sin  /c^x  (C2  sin  ky2  (b  -  y)  +  Z)2  cos  /cy2  (&  -  y))  (C^ e_7n2  +  C“e7n2) 

jweoho 

C2  =  ~ 27n  sin  kxx  (C2  cos  ky2  (b  -  y)  -  D2  sin  ky2  (6  -  y))  (C+e~lnZ  -  C“e7ri2) 
J^oho 

hx 2  =  —  sin  kxx  (C2  sin  ky2  (b  -  y)  +  D2  cos  ky2  (6  -  y))  (C+e_7n2  -  C"e7n2) 
ho 

hy2  =  0 

h^2  =  —  cos  kxx  (C2  sin  ky2  (b  -  y)  +  D2  cos  /cy2  (6  -  y))  (C+e-7"2  +  C“e7712)  (4.21) 
ho 

Boundary  conditions  require  that  the  tangential  electric  field  must  vanish  at  y  =  0 
and  y  =  b,  forcing  the  constants  Ci  and  C2  to  zero.  Continuity  of  tangential  E  and  H 
at  y  =  h  requires  that  kx  and  72n  be  the  same  in  each  subregion  [16].  In  addition,  the 
tangential  components  exl,  ex2,  ezl,  ez2,  exl,  hx2.  hz\  and  hz2  must  be  continuous 
at  the  boundary  between  the  subregions,  y  =  h.  Equating  e2i  and  ez2  of  (4.20)  and 
(4.21),  the  requirement 

ky l'T'n  .  7I’3'  .  i  i  ky2yn  .  7rx  .  .  ±  .  . 

— - sm  —  sm  kvih  ■  B„  =  — - - sin  —  sm  kv2(b  —  h)  ■  G;r  4.22 

jujey  a  3^e0y0  a 
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is  obtained.  Similarly,  the  continuity  between  hz i  and  hz 2  requires 


k  k 

—  cos  (kxx)  cos  (kyih)  ■  B±  =  —  cos  (kxx)  cos  W*  -  A)  '  <7  (4.23) 

/Xl  /io 

Using  either  (4.22)  or  (4.23),  the  propagation  coefficient  of  subregion  1,  B±,  can  be 
expressed  in  terms  of  the  propagation  coefficient  of  subregion  2,  C±,  that  is 


B±  =  -~ky2e1ni  sin  (ky2  (b  -  h))  ±  =  //1  cos  {ky2  (6-  h))  ± 

kyieoHo  sm(kyih)  n  cos  {ky\h) 


(4.24) 


Division  of  (4.22)  by  (4.23)  gives 


sin  kyih  cos  ky2{b  —  h)  =  — —  sin  ky2(b  —  h )  cos  ky\h 
Cl  ^2 


(4.25) 


or,  equivalently 


k  k  9 

—  sin  kyihcos  ky2(b  —  h)  H — —  sin  ky2(b  —  h )  cos  =  0 
6i  62 


(4.26a) 


7m  n  =  +  k2x-kl  =  A&,B  +  k2x-k20  (4.26b) 

such  that  the  wave  numbers  in  each  subregion  satisfy  the  constraint  equation  of 
(4.26b).  Since  both  kyi  and  A72  are  dependent  on  (4.26a)  represents  a  transcen¬ 
dental  eigenvalue  equation  for  possible  values  of  the  PFW  mode-propagation  constant 
7Z.  Therefore,  the  equation  must  be  solved  numerically.  When  the  correct  value  of  7 
is  found,  the  ratio  ^  is  given  by  (4.24).  Incorporating  this  ratio  into  (4.20),  the  field 
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components  of  subregion  1  are 


exi  =  kd'ky2  cos  kxx sin  (C+ e  7nZ  +  Cn  e7"2) 


Gyl  ~ 


jujeofio 

(ko  -  kh) 


sin  kxx  cos  kyiy 


sin  ( kyih ) 
cos  ( ky2  (ft  —  h)) 


cos  (kyih) 


&zl 


jueifio 

-ky lln  ..Sin  (ky 2  (b  -  ft))  , n+ 


(Cie-™  +  C-e-”‘) 


sin  kx  sin  kyly- 


jixeoyo  sin  (ky\h) 

7 n  , _ ..cos (ky2  ( b-h ))  (n+  z 


{pne~lnZ  ~  C~e^) 


hx i  =  —  sin  kxx  cos  kyly 

y0  cos {kyih) 


(C+e-™  -  C-e’-s) 


ky  1  0 

,  .  COS  (/cv2  (ft  -  ft-))  /  „+  „  ,N 

hzi  =  ~  cos  kxx - y— — — - (C+e  7n~  +  Cn  e7n~) 

/io  cos  [kyih) 


(4.27) 


and  the  field  components  of  subregion  2  are 


=  ^xky2  cos sin  (ky2  (ft  -  y))  (<7+ e  7n2  +  Cn  e7"2) 
J^eoho 


ey2 


■  sin  kxx  cos  (ky2  (b  -  y))  (C+ e  7n2  +  Cn  e7nZ) 


jLueoyo 

ez2  =  .  - ‘y27n  sin  kxx  sin  (fty2  (ft  -  2/))  (C+e~7nZ  -  C~elnZ) 

hx 2  =  —  sin  fta-x  cos  (ky2  (ft  -  y))  (C+e_7n2  -  C~e7nZ) 
ho 

fty2  0 


kr 


hz 2  =  —  cos  cos  (ky2  (b  —  y))  (C+  e  7"2  +  Cn  e7"2) 
ho 


(4.28) 


The  wave  number  y2  is  determined  using  the  height  iteration  method  discussed 
in  Section  4.2.6.  The  wave  numbers  ky±  and  ky2,  which  represent  variation  in  the  y 
direction  in  each  subregion,  are  not  expected  to  be  equal  to  each  other.  The  fields  in 
each  of  the  three  regions  having  now  been  constructed,  it  becomes  possible  to  define 
the  system  modes. 


4-2.3  Modal  Analysis.  In  an  empty  waveguide,  only  the  dominant  mode 
propagates  completely  through  the  waveguide.  All  other  higher-order  modes  are 
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z  =  0  z  =  d 


Figure  4.4:  The  TM|0  mode  is  incident  on  the  sample,  and  scatters  the  mode  set 
TM±».  Hybrid  TMy  modes  exist  in  Region  II,  between  2  =  0  and  z  =  d.  In  Region  III 
(z  >  d),  forward  propagating  TMy  modes  exist.  All  modes  for  n  >  0  are  evanescent. 


evanescent  and  rapidly  decay  before  reaching  the  network  analyzer  ports.  For  the 
geometry  described,  which  is  discontinuous  in  y,  the  dominant  mode  is  TMf0  and  the 
scattered  mode  set  is  TM^,  where  the  index  n  indicates  the  mode  number.  Figure 
4.4  shows  a  cross  section  of  the  waveguide,  and  illustrates  the  decay  of  the  evanescent 
modes  below  cutoff  as  well  as  the  propagation  of  the  dominant  mode. 

In  any  region  of  the  system,  the  total  transverse  electric  and  magnetic  field  can 
be  represented  as  the  superposition  of  the  forward  and  reverse  propagating  modes: 


~f^inz 


71=  1 
OO 


H t  =  E  ±an  h„eT7"z 

71=1 


(4.29) 


where  the  coefficients  a±  represent  the  complex  weighting  constants  of  each  mode.  The 
held  vectors  en  and  h„  describe  the  vector  components  of  the  electric  and  magnetic 
fields  of  each  mode,  respectively  [13].  Depending  on  the  value  of  z,  either  the  complex 
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propagation  wave  number  of  free  space  70,  or  that  of  the  PFW  region,  should  be 
used. 

Using  (4.29)  as  a  guide,  the  total  transverse  fields  in  each  region  of  the  system 
can  be  constructed.  In  Region  I  {z  <  0),  a  TM^0  mode  propagates  in  the  forward 
direction,  and  an  infinite  number  of  TMfn  modes  scatter  in  the  reverse  direction  due 
to  scattering  from  the  sample  (see  Figure  4.4).  Now  (4.29)  can  be  expressed  as 

OO 

E  =  a^eioe-70’102  +  E  a^e^e70’1"2 

n=0 

00 

H  =  a^0h10e_7o-lo2:  -  E  a^h^e70’1"*  (4.30) 

71=0 

To2!„  =  (P2  +  Of  )2  -  kl 

Although  an  inhnite  number  of  reverse  scattered  modes  exist,  only  the  n  =  0  mode 
will  propagate  completely  through  the  waveguide  to  the  network  analyzer. 

In  Region  II  (0  <  z  <  d),  the  existence  of  an  inhnite  number  of  forward  and 
reverse  traveling  modes  are  necessary  to  satisfy  the  boundary  conditions  at  the  dis¬ 
continuity  in  the  geometry  of  the  sample.  The  transverse  fields  can  be  written  as 


E 

H 


E  Kn^rne  llnZ 


n= 0 


00 

+  E  Knelne^ 


n= 0 


E  Kn^ine  7l"2  -  E  ^i„hine71™2 


n=0 


n= 0 


(4.31) 


The  use  of  the  tilde  (~)  for  the  Region  II  held  vectors  indicates  their  dependence  on 
position  in  y.  A  different  vector  potential  is  used  to  generate  the  helds  of  each  subre¬ 
gion.  The  propagation  of  modes  in  Region  II  depends  on  several  factors,  including  the 
size  of  the  gap  and  the  nature  of  the  material  under  test.  The  size  of  the  gap  dictates 
whether  a  mode  is  propagating  or  evanescent.  If  the  material  sample  is  lossy,  no  true 
cutoff  exists  and  all  modes  will  propagate,  although  the  modes  that  exist  above  the 
operation  frequency  will  rapidly  decay  [17]. 


In  Region  III  (z  >  d),  the  total  helds  can  be  represented  by  the  forward  prop¬ 
agating  mode  set  of  TMfn;  no  reverse  traveling  waves  are  present  due  to  the  absence 
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of  reflecting  obstacles  in  this  region.  The  final  mode  set  can  be  described  by 


E 

H 


OO 

E  c^nelne"7o’ln^"d) 

n= 0 


E  Clnhlne  d) 


n= 0 


(4.32) 


where  70  is  the  propagation  wave  number  of  free  space.  A  phase  shift  has  been 
introduced  into  the  complex  exponential  to  facilitate  matching  of  boundary  conditions 
between  regions. 


4-2-4  Boundary  Conditions.  It  is  necessary  to  maintain  continuity  of  the 
transverse  E  and  H  held  components  at  the  interfaces  between  each  region,  namely 
at  z  =  0  and  z  =  d,  which  implies  the  satisfaction  of  the  conditions 


4( 

O 

1 

=  ejf  (z  =  0+) 

4( 

z  =  0") 

=  ex  0  =  o+) 

K( 

>  =  0-) 

=  K1  (z  =  0-) 

e11 

y 

{z  =  d~ 

)  =  411  ( z  =  d+ ) 

ei1 

{z  =  d~ 

)  =  elH  (z  =  d+ ) 

K1 

(z  =  d~ 

)  =  h"1  (z  =  d+) 

(4.33) 


Since  the  scattered  mode  set  is  TMy,  the  held  components  hy  do  not  exist.  It  will  be 
shown  that  this  allows  ex  to  be  written  as  a  function  of  the  held  components  ey  and 
hx,  thereby  eliminating  the  need  to  satisfy  ex  explicitly. 

From  (3.1)  and  (3.2),  Maxwell’s  curl  equations  (Faraday’s  Law  and  Ampere’s 
Law)  in  a  source  free  region  are 


V  x  E  =  —juj/j, H 

V  x  H  =  jueE 


39 


Performing  the  curl  in  rectangular  coordinates,  and  letting  hy  =  0,  this  can  be  ex¬ 
panded  into  six  equations, 


dez 

dy 


dey 

dx 


dhx 

dz 


dey  = 
dz 
dex 
dz 

dex  _ 

dy 

dhz 

dy 

_ 

dx 

dhx 

dy 


jujyhx 

(4.34a) 

de 

ir  =  {) 

dx 

(4.34b) 

-jujyhz 

(4.34c) 

=  jueex 

(4.34d) 

=  jujeey 

(4.34e) 

=  jueez 

(4.34f) 

Performing  the  partial  derivative  on  ex  in  (4.34b)  and  using  the  expression  for  ex  from 


(4.17)  yields 


T  ^Iz&x 


dez 

dx 


(4.35) 


where  the  sign  is  chosen  with  respect  to  either  forward  or  reverse  traveling  waves. 
It  will  be  carried  through  the  development.  This  result  is  combined  with  (4.34d)  to 
obtain  the  equality 


1  dez 
=F7 2  dx 


From  (4.34c),  the 


1  dhz 
juje  dy 


(4.36) 


hz 


1  ffey 

—juiy  \  dx 
1  /  dex 

jtvy  V  dy 


_  dex\ 
dy  ) 
dey\ 
dx  ) 


(4.37) 


and  from  (4.34f),  the  expression  for  ez  is 


1  dhx 
jue  dy 


(4.38) 
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Substituting  (4.37)  and  (4.38)  into  (4.36)  and  distributing  the  derivative  operator 
obtains 


1  d  / 

1  dhx\ 

1  d 

\  1 

(  dex 

dey  V 

=F7 zdx  \ 

jue  dy  J 

juje  dy 

K  dy 

dx  J 

±1  d2hx  _  1  /  d2ex  _  d2ey  \ 

juje'fz  dxdy  co2ey  \  dy2  dydx  J 


(4.39) 


Taking  the  two  terms  on  the  right  hand  side  of  (4.39),  substituting  =  —k2ex, 
multiplying  both  sides  by  u;2e/i,  and  rearranging  terms  yields 


1  /  ±ju>n  d2hx  d2ey  \ 
k2  \  dxdy  dydx ) 


(4.40) 


The  total  transverse  field  consists  of  both  forward  and  reverse  propagating  modes, 
and  the  (±)  is  applied  to  each  type  of  mode  respectively.  At  material  boundaries, 
regardless  of  the  presence  of  discontinuities,  tangential  E  and  H  and  their  derivatives 
must  be  continuous  [10,16].  Having  shown  that  ex  is  linearly  dependent  on  the  second 
derivatives  of  ey  and  hx,  the  interfacial  boundary  conditions  on  ey  and  hx  will  be 
satisfied  analytically.  The  resultant  continuity  of  ex  can  be  shown  numerically  given 
a  successful  solution  of  the  mode  matching  matrix  of  (4.56). 

Region  11  is  divided  into  two  subregions,  each  with  different  material  parame¬ 
ters.  Therefore,  continuity  of  ey  and  hx ,  and  also  (4.40),  must  satisfy  the  boundary 
conditions  (4.33)  piecewise,  that  is 

e7| 

y 1 2=0- 

e1 1 

Hz=o- 


ey\  ■, 0  <y  <  h 
e'yiX  <V  <b 


(4.41) 


z= 0+ 


ex\ ,  0  <  y  <  h 
eTxi,h  <  y  <b 


(4.42) 


z=0+ 
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h1 1 

x\z=0 


K[,0  <y  <h 
hx2,h  <  y  < b 


(4.43) 


2=0  + 


at  the  z  =  0  plane.  Additionally,  at  the  z  =  d  plane,  the  piecewise  equalities 

(4.44) 


ey{l  °  <y  <h 


eTyi,h  <  y  <  b 

4i,o  <  y  <  h 
eiih  <  y  < b 

hii,  o  <  y  <  h 
hx2,h  <  v  < b 


=  e ; 


mi 

y  I  z=d+ 


z=d~ 


=  e 


mi 


z=d+ 


z=d~ 


=  h 


mi 

31  \z=d+ 


(4.45) 


(4.46) 


z=d~ 


must  hold,  where  e/J ,  e77  and  h77  are  functions  of  position  in  y,  corresponding  to  the 
vector  field  components  of  each  respective  subregion.  The  total  fields  are  matched 
at  the  boundary  are  the  superposition  of  an  infinite  number  of  higher  order  modes, 
which  must  be  truncated  through  practically  to  N  modes  for  computational  purposes. 

Now  that  the  required  transverse  field  components,  namely  ey  and  hx,  have 
been  identified,  the  boundary  conditions  of  (4.33)  can  be  applied  to  the  mode  sets  of 
(4.30),  (4.31),  and  (4.32).  The  superposition  of  N  modes  yields  the  transverse  fields 
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(4.47) 


n= 0 


n= 0 


n= 0 
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where  the  leading  coefficients  have  been  normalized  with  respect  to  af0,  resulting  in 


aln 
-f-  ^ 

r  —  hn 
n  „+  ’ 

,  b+n 

ln  „+  ’ 

r+ 

rr ?  1  In 

±Tl  —+ 

(4.48) 

a10 

a10 

a10 

a10 

The  constants  of  (4.48)  correspond  to  the  interfacial  reflection  and  transmission  co¬ 
efficients  of  the  sample,  at  the  front  (Tn,tn)  and  back  ( rn,Tn )  sample  interfaces, 
respectively.  The  required  S-parameters  of  the  system  can  be  determined  from  the 
mode  coefficients,  as  they  are  simply 


(4.49) 


If  the  sample  is  a  simple  material,  the  reverse  S-parameters  can  be  accurately  equated 
to  the  forward  parameters,  that  is 


gthy 

r-thy 
°21 

4-2.5  Making  a  Well-Posed  Problem.  The  system  of  (4.47)  contains  4  equa¬ 
tions  and  4 N  unknowns  and  is  therefore  underdetermined  for  n  >  1.  If  the  problem  is 
to  be  well-posed,  an  equal  number  of  equations  and  unknowns  must  exist.  Since  the 
number  of  unknowns  cannot  be  reduced,  a  method  to  construct  AN  equations  must 
be  used.  Physically,  this  matrix  represents  the  coupling  of  modes  between  regions 
of  empty  waveguide  and  partially  filled  waveguide.  While  only  the  dominant  mode 
( n  =  0)  will  propagate  through  the  entire  guide,  complex  power  will  also  be  transmit¬ 
ted  and/or  stored  by  the  higher-order  modes  (n  >  0)  which  must  be  accounted  for  if 
an  accurate  model  of  the  fields  is  to  be  proposed. 

The  satisfaction  of  the  boundary  conditions  between  regions  presented  in  Section 
4.2.4  required  equating  a  single  vector  component  of  both  the  transverse  electric  and 
magnetic  fields,  that  is  the  ey  and  hx  components.  A  suitable  testing  operator  should 


(4.50) 
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reduce  system  complexity  without  altering  the  resultant  S-parameters.  To  ease  the 
required  computations,  mutually  orthogonal  operators  are  desirable.  Inspection  of 
(4.17),  (4.20)  and  (4.21)  shows  that  the  sin  (kxx)  term  is  common  to  the  required 
vectors.  The  argument  of  the  cosine  term,  however,  is  different  in  each  region.  Since 
it  is  known  that  in  Region  I,  ky  =  an  appropriate  testing  operator  is 

0  =  cos  P  —  0, 1,  2,  ...N  (4-51) 

to  be  applied  to  each  mode  vector, 


(0,  f)  =  /  fp  ■  fdy  (4.52) 


which  is  the  inner  product  of  a  mode  vector  /  and  the  testing  function  0.  This 
choice  of  operation,  when  used  on  the  fields  of  Regions  I  and  III,  will  cause  many 
terms  to  vanish,  since  sinusoidal  functions  of  different  mode  indices  are  mutually 
orthogonal.  The  linear  dependence  of  ex  on  the  field  components  ey  and  hx  removes 
the  requirement  to  integrate  the  testing  function  over  the  waveguide  cross  section,  as 
is  often  done  in  other  methods  [7,13,17,18].  The  testing  operation  is  applied  term  by 
term  to  the  linear  system  of  equations  in  (4.47),  which  can  then  be  written  as 

N- 1  N- 1  N- 1 

Y  Tn  (eln,  0)  -  Y  Tn  ^ln’  “  Y  tn  ^ln’  ^)  +  0  =  (  e10)  0) 

71=0  71=0  71=0 

TV— 1  N—l  N—l 

Y  rw  (hlrt,  0)  -  Yrn(hln,^  +  Y.tn  (hln,0),+O  =  (hlO,0) 

71=0  71=0  71=0 

N  N-l  N—l 

0  +  Y Vn  ^ ellnd + Y tn  ^in’  ^  e~llnd  _  Y  Tn  ^ein)  ^  =0 

71=0  71=0  71=0 

N  N-l  N-l 

0  +  J>  ^hin,  0^>  ellnd  -  y  tn  (hin,  0^  e~llnd  +  Y  Tn  (hi„,  0)  =0  (4.53) 

71=0  71=0  71=0 

To  further  reduce  the  required  number  of  calculations,  it  is  useful  to  observe  that 
the  transverse  e  and  h  field  components  have  similar  spatial  variation,  containing 
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harmonic  functions  of  identical  arguments.  This  means  that  the  testing  operation 
integrals  need  only  be  calculated  once,  and  then  can  be  multiplied  by  an  appropriate 
scaling  factor  to  correspond  to  the  necessary  field  component.  To  facilitate  anticipated 
matrix  algebra,  the  testing  operations  of  (4.53)  corresponding  to  the  electric  field 
mode  vectors  can  be  represented  by  submatrices,  such  that 


Mnp 

Nnp  (,  •  Unp  T  Vnp 
Unp  (ei ,  'ifj ) 

Uip  (®2j 


(4.54) 


where  Q  is  the  scaling  factor  from  (4.24).  Submatrices  Unp  and  Vnp  correspond  to  the 
testing  operation  applied  to  the  electric  held  vector  of  subregions  1  and  2,  respectively. 
When  the  inner  product  is  performed  on  the  magnetic  held  vectors,  the  resultant 
submatrices  are 

Pnp 

Qnp  =  C  ‘  Pnp  T  Gnp 

Pnp 

Gnp  =  (e2,  V’) 

where  Z  is  the  ^-directed  wave  impedance  in  each  region,  and  the  notation  is  that  of 
(4.54).  The  respective  integrals  of  (4.54)  and  (4.55)  are  presented  in  Appendix  B 


Using  (4.54)  and  (4.55),  the  system  of  (4.53)  can  be  represented  in  matrix  form 


as 


Mnv  —Nnp 

~  Pnp  Qnp 

0  e~llntNnp 

0  e-^Qnp 


- Nnp 

0 

Qnp 

0 

pUnt  AT 
°  1  ynp 

- Mnp 

—  pllntr) 

c  ^5 np 

—  P 

1  np 

r„ 

i 

_ 1 

rn 

—  Pop 

^n 

0 

T 

- L  n 

0 

(4.56) 


The  system  is  now  in  the  form  Ax  =  B,  where  each  submatrix  of  A  is  of  dimension 
N  x  A,  making  the  entire  matrix  4 N  x  AN.  Both  the  solution  set  x  and  the  excitation 
vector  B  are  4 N  x  1  column  vectors. 
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In  order  to  numerically  verify  the  continuity  of  the  ex  fields  across  the  interfa¬ 
cial  boundary,  the  coefficients  of  the  solution  set  x  are  used  to  check  the  boundary 
conditions  on  ex, 


N-l  /  N—l  N—l  \ 

^ n&x,ln  (  Y'-  tnCx  \n  T  'Y.  J  0 

\ra=0  n= 0  / 

N-l  \  N-l 


n= 0 
/  N—l 

Y  infix, ine~lzd  +  J]  rnex,inelzd  )  -  J]  TneXj ln  =  0 

\  n=0  n=0  J  n=0 


(4.57) 


This  was  calculated  with  several  combinations  of  height  h  and  total  modes  N,  with 
absolute  differences  on  the  order  of  ICR15.  Together  with  (4.40),  this  verifies  that  the 
fields  are  continuous  across  the  regional  boundaries. 

As  has  already  been  indicated,  the  S-parameters  S^y  and  S^y  are  identically 
Tj  and  7\  of  the  solution  set.  The  minimization  equation  (4.1)  can  now  be  used  to 
search  for  the  correct  values  of  e  and  /i. 


4-2.6  Height  Iteration  Method.  It  has  been  stated  that  the  eigenvalue  equa¬ 
tion  of  (4.26a)  is  a  transcendental  equation  for  possible  values  of  the  complex  propa¬ 
gation  wave  number  ^z^pfw  and  must  be  solved  numerically  [1].  A  good  initial  guess 
is  imperative  to  successfully  finding  the  root  of  the  equation.  The  value  of  7 z,ffw  of 
the  fully  hllwed  waveguide  (FFW)  is  used  as  the  initial,  unperturbed  guess  for  7 z,pfw 
in  (4.26a)  with  an  initial  guess  as  h  =  b  —  5,  where  5  is  a  small  value. 

A  Newton-Raphson  root  search  algorithm,  using  central-difference  derivatives, 
uses  these  values  to  converge  on  a  new  7,  which  becomes  the  initial  guess  in  the  next 
iteration.  The  height  h  is  decreased  again  by  5,  and  the  process  repeats  until  the 
input  height  h  is  the  actual  value.  If  the  final  iteration  converges,  then  7  is  accepted 
as  7 pfwi  which  is  used  to  fill  the  mode  matrix  of  (4.56).  The  top-level  root  search, 
using  the  minimization  equation  of  (4.1),  continues  the  iteration  for  parameters  e  and 
A 


46 


4- 3  Reference  Plane  Independent  Measurement 

The  two  measurements  discussed  in  the  previous  section,  5)  i  and  S2i,  are  the 
reflection  and  transmission  S-parameters,  respectively.  The  reflection  measurement  is 
highly  dependent  on  the  sample  location  relative  to  the  calibration  plane  along  the 
z-axis  (see  Figure  4.1),  since  a  change  in  position  incurs  a  two-way  phase  delay  (or 
advance)  in  the  measurement.  The  transmission  measurement  is  not  sensitive  to  the 
placement  of  the  sample. 

To  perform  calculations,  it  is  necessary  to  extract  the  true  S-parameters  (de¬ 
noted  with  superscript  s)  from  the  measured  S-parameters  (denoted  with  superscript 
ms )  obtained  from  the  network  analyzer.  This  section  will  assume  that  the  sample 
has  thickness  d,  and  is  in  a  waveguide  sample  holder  of  width  w.  Also,  the  notation 
kz,  instead  of  jz,  is  used  for  the  propagation  wave  number.  Using  Figure  4.5  (a)  as  a 
guide,  and  using  complex  exponential  notation  to  denote  phase  shifts,  the  measured 

5- parameters  are: 


Qms  _  QS 

°11  —  °11 

s™s  =  e-Mu,-d)  s^e-jk(w-d)  _  e-jk2{w-d) £|2 
=  e 

Cjms  _  g—jk(w—d)  Cjs 


(4.58) 


Upon  de-embedding  the  S-parameters,  the  actual  material  S-parameters  are  found  to 
be 


as 

°22 

as 

°2l 

as 

°12 


Cms 

°11 

ejk2(w-d)Sms 
ejk(w-d)  Qjms 
ejk(w-d)  Cjms 


(4.59) 


The  preceding  equations  presume  that  the  sample  is  precisely  aligned  with  the 
calibration  plane.  If,  however,  the  alignment  is  erroneous,  or  if  the  sample  shifts  during 
actual  measurements,  the  extraction  of  parameters  in  (4.59)  is  invalid.  Consider  the 
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(a) 


Figure  4.5:  Diagram  of  a  waveguide  sample  undergoing  a  longitudinal  position  shift, 
(a)  Sample  (s)  and  measured  (ms)  S-parameters  of  a  sample.  The  sample  is  aligned 
with  Calibration  Plane  1.  (b)  The  sample  has  been  shifted  5  along  the  axis,  and 
ii  <  i.  The  two  reflection  measurements  now  contain  an  additional  phase  shift  of 

e±jfc2<5 
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situation  of  Figure  4.5  (b)  when  the  sample  is  shifted  along  the  waveguide  a  small 
distance  S  <C  d.  The  actual  S-parameters  are  then 

gs  gjk25  Qms 

_  gjk2(w—d—5)cjms 

5|j  =  ejk(w~d)  SZl3 

Qs  _  gjk(w—d)  Qms 

It  can  be  seen  that  the  forward  and  reverse  transmission  measurements  of  (4.60) 
and  (4.59)  are  the  same.  However,  the  forward  and  reverse  reflection  measurements 
of  (4.60)  incorporate  the  two-way  phase  delay/advance  caused  by  the  shift.  Even  a 
small  value  of  <5  will  cause  a  significant  phase  shift  of  the  S-parameters. 

It  is  quite  possible  that,  in  the  course  of  handling  the  material  sample  and  the 
waveguide  system,  the  exact  distance  between  the  sample  and  the  calibration  plane 
may  not  be  known.  During  high-temperature  measurements,  this  problem  can  be 
exacerbated  by  waveguide  expansion  in  extreme  heat.  To  perform  accurate  material 
characterization  in  these  situations,  it  is  necessary  to  eliminate  the  dependence  on 
the  reference  plane.  By  multiplying  the  forward  and  reverse  reflection  measurements 

11^22  =  (e  Wl  )  \e  b22  ) 

„jk2(w—d)  oms  nms 

=  eJ  Jbn  <b22 

and  the  forward  and  reverse  transmission  measurements 

ns  os  ( ^jk(w-d)  oTns\  ( ^jk(w-d)  oms\ 

^21^12  =  (e  ^21  )  { e  ) 

=  ej2k{w~d)  S™  S™  (4.62) 

the  dependence  on  5  drops  out.  The  remaining  quantities  (w,d,  and  kz)  are  all  known 
values.  When  using  a  root  search  method  to  determine  permittivity  and  permeability, 


(4.61) 


(4.60) 
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these  products  can  be  used  in  a  system  of  minimization  equations 


rithy  rithy  riexp  riexp 

‘-’ll'  ‘-,22  —  ‘-’ll  D22 

rithy  rithy  riexp  riexp 
‘-’21'  ‘-’12  —  ‘-’21  ‘-’12 


<  tol 

<  tol 


(4.63) 


In  a  simple  material,  the  respective  forward  and  reverse  reflection  and  transmission 
measurements  are  equal,  i.e.  Sn  =  S22  and  S21  =  S\2-  This  allows  Id  and  Tf  of  the 
solution  set  of  (4.56)  to  be  used  for  both  the  forward  and  reverse  coefficients,  saving 
significant  computational  effort. 


4-4  Summary 


A  method  for  characterizing  an  electromagnetic  material  that  partially  fills  a 
rectangular  metal  waveguide  in  one  dimension  was  presented  using  TMy  modal  anal¬ 
ysis.  The  matrix  A  contains  the  region-to-region  mode  coupling  information  and  the 
solution  vector  x  relates  A  to  the  held  excitation  vector  B.  The  theoretical  reflec¬ 
tion  and  transmission  coefficients  of  the  dominant  mode,  S’j1(y  and  S^y ,  are  extracted 
from  x  and  compared  to  the  experimental  S-parameters.  A  2-D  Newton-Raphson 
root  search  iterates  the  parameters  er  and  yr  until  the  absolute  difference  between 
the  theoretical  data  and  experimental  data  is  within  a  specified  tolerance. 


In  addition,  a  method  was  developed  for  removing  measurement  dependence  on 
a  sample’s  axial  position  in  the  waveguide  with  respect  to  the  reference  plane.  This 
was  accomplished  by  multiplying  the  forward  and  reverse  reflection  measurements 
^S,)111>  and  and  forward  and  reverse  transmission  measurements  f  and  S^y 
which  are  then  compared  to  the  products  of  the  respective  experimental  S-parameters. 

A  2-D  Newton-Raphson  root  search  is  again  used  to  arrive  at  a  usable  solution. 
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V.  Results 


The  PFW  measurement  correction  and  the  reference  plane  independence  analysis 
were  developed  in  the  previous  chapter  for  the  case  of  a  single  air  gap  between  the 
top  of  the  material  sample  and  the  waveguide.  The  analysis  uses  mode  matching  to 
calculate  the  full  set  of  theoretical  scattering  parameters  and  compares  them  to  the 
experimentally  measured  scattering  parameters.  A  two-dimensional  Newton-Raphson 
root  search  is  used  to  minimize  the  difference  between  the  theoretical  and  experimental 
S-parameters  using  the  equations 

S'n  (w,  e,  n)  -  Snp  (u) 

^2iy  -  ^2iP  M 

where  it  is  accepted  that  the  actual  values  of  permittivity  e  and  permeability  /i  will 
drive  the  two  functions  of  (5.1)  below  the  specified  tolerance. 

5.1  Test  Procedure 

Microwave  measurements  were  performed  at  room  temperature  in  two  frequency 
bands,  S-band  (2.6  —  3.95  GHz)  and  X-band  (8.2  —  12.4  GHz).  The  room  tempera¬ 
ture  characterizations  were  done  on  samples  with  deliberately  machined  gaps  in  order 
to  simulate  the  effect  of  a  high  temperature  waveguide  expansion.  Measurement  of  the 
physical  dimensions  of  the  samples  was  performed  in  conjunction  with  the  microwave 
measurements. 

Room  temperature  tests  of  samples  machined  with  a  gap  sufficiently  mimicked 
the  conditions  of  a  high  temperature  test  to  the  satisfaction  of  the  research  sponsor. 
All  data,  both  at  S-  and  X-band,  was  collected  on  a  Hewlett-Packard  8510C  Network 
Analyzer  located  at  the  AF1T  Microwave  Lab  between  the  months  of  November  2006 
and  February  2007.  The  measurement  apparatus  for  S-  and  X-band  can  be  seen 
in  Figure  5.1  and  Figure  5.2,  respectively.  A  Thru-Reflect-Line  (TRL)  calibration 
scheme  was  used  to  calibrate  out  the  sample  holder. 


<  tol 

<  tol 


(5.1) 


51 


The  results  presented  are  the  characterization  of  two  materials,  acrylic  and 
FGM-125.  The  acrylic  material  is  a  commercially  available,  ideally  lossless  dielectric. 
FGM-125  is  a  commercially  available  rubberized  magnetic  shielding  material,  having 
both  electric  and  magnetic  losses,  manufactured  by  Emerson  &  Cuming.  Measure¬ 
ments  of  the  acrylic  were  only  done  at  S-band,  while  FGM-125  was  measured  at 
S-band  and  X-band.  This  research  is  intended  specifically  for  magnetic  materials,  but 
it  can  be  applied  generally,  hence  the  inclusion  of  the  acrylic  measurements. 

Parameter  data  is  presented  visually  using  the  convention  er  =  e'  —  je"  and 
Hr  —  n'  —  jh".  Truth  data  was  calculated  using  NRW  on  a  measurement  of  a  sample 
that  completely  Elled  the  waveguide.  RPI  was  not  used  to  calculate  the  truth  data, 
since  it  is  not  a  standard  measurement  technique.  The  imcorrected  data  is  the  raw 
data  of  the  PFW  measurement,  and  NRW  is  used  to  extract  the  effective  permittivity 
and  permeability  of  this  measurement.  Additionally,  the  minimization  equations  of 
(4.63)  are  used  in  conjunction  with  PFW  theory  to  characterize  samples  using  the 
reference  plane  independent  formulation.  The  corrected  data  corresponding  to  this 
method  is  labeled  RPI.  The  RPI  correction  is  always  calculated  with  the  same  number 
of  modes  as  the  standard  mode-matching  correction. 

5.1.1  Error  Analysis.  Throughout  the  course  of  the  research,  every  attempt 
was  made  to  minimize  conceivable  sources  of  error.  The  well-known  TRL  technique 
was  used  to  calibrate  the  waveguide/HP  8510  system,  minimizing  systematic  instru¬ 
mentation  uncertainty.  The  waveguides  were  mounted  on  fixed  racks  to  minimize 
unnecessary  movement  of  cables  and  prevent  accidental  shifting  of  the  sample  in  the 
sample  holder  (see  Figures  5.1  and  5.2).  The  use  of  precision  alignment  pins  between 
the  waveguide  flanges  and  sample  holder  ensured  consistency  of  the  system  during 
calibration  and  measurements. 

One  unavoidable  source  of  error  was  the  measurement  of  the  physical  dimen¬ 
sions  of  each  sample  under  test.  Using  calipers  with  ±2  mil  accuracy,  each  sample 
dimension  (width,  height,  thickness)  was  measured  five  times,  and  the  average  value 
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Figure  5.1:  Mounted  S-band  waveguide  attached  to  HP  8510  Network  Analyzer. 
The  mounting  rack  minimizes  unnecessary  movement  of  the  cables. 


Figure  5.2:  Mounted  X-band  waveguide  attached  to  HP  8510  Network  Analyzer 
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was  recorded.  A  differential  error  analysis  was  performed  to  determine  the  effect  of 
uncertainty  in  height  h  and  thickness  d  on  the  final  calculation  of  permittivity  and 
permeability. 

It  has  been  suggested  that  the  real  and  imaginary  components  of  permittivity 
and  permeability  be  described  as  functions  of  h  and  d  [12,13].  However,  this  choice  of 
functions  does  not  reflect  the  mechanics  of  the  algorithm  itself,  since  the  actual  solu¬ 
tion  set  is  theoretical  S-parameters.  An  error  analysis  of  the  S-parameters  themselves 
would  not  be  particularly  useful,  since  the  purpose  of  the  research  is  to  calculate  e  and 
(jl.  Therefore,  the  error  analysis  of  e  and  /i  will  consider  propagation  of  uncertainty 
in  the  calculated  S-parameters  due  to  the  measurement  uncertainty  in  sample  height 
and  thickness. 

Let  k  be  either  relative  permittivity  or  permeability,  such  that 

k  —  k'  (h,  d)  —  jn"  (h,  d) 


The  real  part  of  the  calculated  material  parameter  due  to  a  height  h  and  thick¬ 
ness  d  is  expanded  in  a  Taylor  series  around  the  point  h0  and  d0, 


//,  ,  n  ,  ,  r «  _  ,/,  7  \  .  dK'(h0,d0)  8k1  (ho,  do) 

k  (ho  -\-  oh,  do  T  od)  —  k  (ho,  do)  T - — - oh  T - — ; - od  T 


dh 


dd 


(5.2) 


where  5h  and  5d  are  the  measurement  uncertainties  (±2  mils)  in  height  and  thickness, 
respectively.  The  higher  order  terms  of  the  Taylor  series  can  be  neglected  since  the 
uncertainty  is  so  small.  The  total  uncertainty  in  k'  can  now  be  approximated  by 


<W  =  k!  (ho  +  5h,  do  +  5d )  —  n'  (h0,  do) 


8k'  (h0,  d0)  xu  8k' (h0,  do) 


8h 


-5h 


dd 


5d 


(5.3) 


which  represents  the  uncertainty  contribution  of  both  measurements.  Since  an  ana¬ 
lytic  expression  for  the  partial  derivatives  does  not  exist,  they  must  be  calculated  via 
a  numerical  approximation.  Using  forward  differences,  the  two  partial  derivatives  of 
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(5.3)  are 


3k!  {hp,  dp) 
dh 


k'  (ho  +  Sh,  do) 


Sh 


k'  (h0ld0) 


dd  (ho,  do) 

dh  6h 


k'  (ho,  do  +  Sd) 


Sd 


The  quantities  Sh  and  Sd  can  in  general  be  positive  or  negative,  but  this  analysis 
will  proceed  considering  the  compounded  uncertainty  of  overmeasuring  both  sample 
dimensions.  While  all  the  errors,  assuming  they  are  independent  and  random,  could 
be  combined  in  quadrature,  a  worst-case  approximation  is 


|fa'l  =  1*4  +  *41  <  1*41  + 1*4 


(5.4) 


The  Triangle  Inequality  shows  that  the  absolute  value  of  the  sum  of  the  errors  is  always 
less  than  the  sum  of  the  absolute  values  [25].  Substituting  the  real  and  imaginary 
parts  of  e  and  /i  in  for  k  completes  the  error  analysis. 

The  height  iteration  method  used  to  solve  the  eigenvalue  equation  for  the  correct 
value  of  the  complex  propagation  constant  in  the  PFW  region  is,  unfortunately, 
not  consistently  stable.  Providing  an  initial  guess  of  7 ffw  for  a  particular  mode  does 
not  guarantee  that  the  root  search  for  7 pp\y  will  converge  on  the  next  consecutive 
root,  if  it  converges  at  all.  It  is  also  possible  that  the  root  search  will  find  the  same 
root  twice,  or  skip  a  root.  If  this  happens,  the  subsequent  filling  of  the  mode  matrix 
in  (4.56)  will  be  inaccurate,  due  to  the  use  of  incorrect  values  of  Verification  of 
sample  homogeneity  is  also  beyond  the  scope  of  this  research,  so  it  is  assumed  that 
material  samples  are  homogeneous. 

5.2  Room  Temperature  PFW  Results 

5.2.1  Acrylic.  Although  this  mode-matching  technique  is  meant  to  be  ap¬ 
plied  to  magnetic  materials,  a  correction  can  be  applied  to  a  dielectric  material  in 
order  to  verify  that  the  algorithm  is  working.  Therefore,  acrylic  samples  were  ob- 
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tained  from  the  AFIT  machine  shop,  manufactured  to  the  dimensions  shown  in  Table 
5.1.  Acrylic  is  non-magnetic,  so  it  was  assumed  that  the  sample  permeability  is  that 


Table  5.1:  Acrylic  Samples,  S-Band 


Sample  # 

Thickness  (in) 

Width  (in) 

Height  (in) 

Effective  Gap  (  ±2  mils) 

1 

0.24 

2.84 

1.340 

0 

2 

0.24 

2.84 

1.328 

12 

of  free  space.  The  PFW  correction  was  performed  on  Sample  2,  measured  to  have 
an  effective  air  gap  of  12  mils.  A  correction  using  5  modes  was  performed,  with 
the  results  for  real  permittivity  displayed  in  Figure  5.3.  Acrylic  is  essentially  lossless 
(e"  ~  0),  so  the  imaginary  component  is  not  presented  graphically.  The  corrected  per¬ 
mittivity  is  within  2%  of  the  true  value  for  the  entire  band,  and  provides  a  noticeable 
improvement  over  the  uncorrected  data.  When  combined  with  the  RPI  formulation 
(also  using  5modes),  the  correction  is  nearly  perfect. 

Having  shown  the  performance  of  the  correction  for  this  simple  case,  it  is  now 
appropriate  to  consider  the  more  complex  situation  of  magnetically  lossy  material. 

5.2.2  FGM-125.  Samples  of  FGM-125  were  obtained  and  machined  to  the 
dimensions  shown  in  Table  5.2,  to  be  used  in  tests  at  S-band  (2.6  -  3.95  GHz).  In 
addition,  samples  were  machined  to  the  dimensions  shown  in  Table  5.3  to  be  used  in 
tests  at  X-band  (8.2  -  12.4  GHz).  The  S-band  results  will  be  discussed  first. 


Table  5.2:  FGM-125  Samples,  S-Band 


Sample  # 

Thickness  (±2  mils) 

Width  (in) 

Height  (in) 

Effective  Gap  (±2  mils) 

1 

125 

2.84 

1.340 

0 

2 

125 

2.84 

1.313 

27 

3 

125 

2.84 

1.259 

81 

Sample  1,  which  completely  filled  the  waveguide  cross  section,  was  used  as  the 
reference  sample.  A  10  mode  correction  was  applied  to  Sample  2,  measured  with  a 
27  mil  top  air  gap.  In  Figures  5.4  and  5.5,  it  can  be  seen  that  while  the  true  complex 
permittivity  is  reasonably  constant  across  the  frequency  band,  permeability  is  not. 
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- 5  Modes - RPI  - True  (NRW) . Uncorrected 


Figure  5.3:  Measured  real  permittivity  of  acrylic  with  12  mil  air  gap.  Acrylic  is  both 
lossless  and  non- magnetic,  so  it  is  assumed  e"  =  0  and  fir  —  1.  A  5  mode  correction 
provides  a  significant  improvement.  The  addition  of  the  RPI-technique  results  in 
a  nearly  perfect  correction.  Error  bars  indicate  ±2  mil  uncertainty  in  height  and 
thickness. 
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Table  5.3:  FGM-125  Samples,  X-Band 


Sample  # 

Thickness  (±2  mils) 

Width  (in) 

Height  (in) 

Effective  Gap  (±2  mils) 

1 

125 

0.90 

0.40 

0 

2 

125 

0.90 

0.39 

10 

3 

125 

0.90 

0.35 

45 

This  is  an  expected  property  of  magnetic  materials,  that  both  p!  and  //'  tend  to 
decrease  quickly  for  increasingly  high  frequencies  [21]. 

From  Figure  5.4  (a),  it  is  apparent  that  the  e'  extraction  is  highly  sensitive  to 
gap  size,  since  the  uncorrected  measurement  is  consistently  8%  below  the  true  data. 
The  modal  correction,  both  with  and  without  the  RP1  formulation,  converges  to  a 
slightly  higher  value  of  real  permittivity,  just  outside  the  range  of  the  error  bars, 
within  5%  of  the  truth  data.  The  uncorrected  data  for  e"  is  sufficiently  close  to  the 
true  value  of  zero,  so  a  modal  correction  is  not  actually  necessary.  However,  it  can  be 
seen  that  inclusion  of  RP1  yields  a  superior  result  than  when  it  is  neglected. 

Examining  the  plots  of  real  and  imaginary  permeability  in  Figure  5.5  it  can  be 
seen  that  these  measurements  are  not  as  sensitive  to  top/bottom  air  gap,  given  the 
proximity  of  the  uncorrected  data  to  the  truth  data.  The  10  mode  correction  does 
not  much  improve  the  extraction  for  either  real  or  imaginary  parts,  although  using 
RPI  has  a  considerable  effect,  especially  for  the  imaginary  component. 

In  a  lab  environment,  it  is  expected  that  the  waveguide  samples  will  not  be 
deliberately  machined  with  large  gaps.  Even  in  a  high  temperature  measurement, 
unless  the  sample  falls  over  onto  its  side,  extremely  large  gaps  are  not  expected. 
However,  in  order  to  illustrate  a  more  drastic  case  of  a  PFW,  and  to  test  the  accuracy 
of  the  modal  method  for  large  gaps,  Sample  2,  with  a  gap  of  81  mils  was  tested.  As  the 
largest  gap  size  under  test,  it  will  be  easier  to  see  the  benefit  of  including  successive 
higher-order  in  the  mode-matching  correction.  Therefore,  a  10  mode  correction  is 
presented  in  Figure  5.6  and  Figure  5.7. 
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- 10  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


- 10  Modes - RPI  - True  (NRW) . Uncorrected 


Figure  5.4:  Permittivity  of  FGM-125  in  a  PFW  with  27  mil  top  air  gap.  A  10  mode 
solution  overcorrects  the  (a)  real  part,  and  yields  an  exact  match  between  3  and  3.1 
GHz  for  the  (b)  imaginary  part.  The  inclusion  of  RPI  is  within  the  error  bounds  in 
(a),  but  gives  superior  performance  in  (b).  Error  bars  indicate  ±2  mil  uncertainty  in 
height  and  thickness. 
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- 10  Modes  -  -  -RPI  - True  (NRW) . Uncorrected 


- 10  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


Figure  5.5:  Permeability  of  FGM-125  in  a  PFW  with  27  mil  top  air  gap.  Correction 
using  10  modes  has  almost  no  difference  from  the  uncorrected  data  for  either  the  (a) 
real  or  (b)  imaginary  components.  The  RPI,  however,  is  an  improvement  to  both 
components,  and  is  nearly  exact  for  //'.  Error  bars  indicate  ±2  mil  uncertainty  in 
height  and  thickness. 
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- 10  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


- 10  Modes - RPI  - True  (NRW) . Uncorrected 


GO 

Figure  5.6:  Permittivity  of  FGM-125  in  a  PFW  with  81  mil  top  air  gap.  10  mode 
correction  applied,  (a)  Use  of  RPI  is  within  the  error  bounds  of  the  10  mode  cor¬ 
rection  for  e' .  (b)  All  extractions  are  reasonably  similar.  Error  bars  indicate  ±2  mil 
uncertainty  in  height  and  thickness. 
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- 10  Modes  -  -  -RPI  - True  (NRW) . Uncorrected 


- 10  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


Figure  5.7:  Permeability  of  FGM-125  in  a  PFW  with  81  mil  top  air  gap.  10  mode 
correction  applied,  (a)  Neither  correction  improves  upon  the  raw  data.  (b)Use  of  RPI 
gives  better  results  than  the  modal  correction  alone,  nearly  perfect  alignment  with 
the  truth  data.  Error  bars  indicate  ±2  mil  uncertainty  in  height  and  thickness. 
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The  correction  of  e'  in  Figure  5.6  (a)  yields  excellent  results,  as  the  truth  data 
is  within  the  uncertainty  of  the  modal  correction,  both  with  and  without  the  RPI 
formulation.  Considering  the  correction  to  e"  in  Figure  5.6  (b),  once  again  the  uncor¬ 
rected  data  yields  an  acceptable  value.  Of  the  modal  correction  methods,  however, 
RPI  is  the  preferred  solution  ,  since  it  provides  the  most  improvement.  Likewise,  in 
Figure  5.7,  the  use  of  RPI  is  preferred  as  a  correction  to  //,  and  y!' . 

The  difference  between  the  uncorrected  permittivity  data  of  the  27  mil  gap  (Fig¬ 
ure  5.4  (a))  and  the  81  mil  gap,  and  the  relative  stability  of  the  other  three  complex 
parameters,  is  a  further  indication  of  the  sensitivity  of  permittivity  extractions  to  this 
particular  PFW  geometry.  This  is  due  to  the  field  pattern  in  the  PFW  region. 

The  electric  field  vector  ey  is  discontinuous  in  y  between  the  material  and  free 
space,  creating  a  capacitive  charge  distribution  [22].  The  measurement  of  electric 
permittivity,  therefore,  is  very  sensitive  to  this  geometry.  However,  the  magnetic  field 
vector  hx,  for  a  given  value  of  ?/,  is  continuous  in  x.  Accordingly,  the  measurement 
of  magnetic  permeability  is  much  more  stable  than  the  permittivity  measurement  for 
top  air  gaps.  If,  for  example,  the  air  gap  was  located  in  the  other  dimension  (i.e. 
left/right),  an  inductive  charge  distribution  would  be  created.  It  can  reasonably  be 
assumed  that  in  this  scenario,  the  permeability  measurement  would  suffer  more  than 
permittivity,  using  similar  reasoning. 

The  second  round  of  room  temperature  measurements  were  performed  at  X- 
band,  using  the  samples  described  by  Table  5.3.  The  parameters  extracted  using 
NRW  on  measurements  of  Sample  1,  the  fully  filled  case,  are  assumed  to  be  the  truth 
data.  The  first  measurements  were  taken  of  Sample  2,  with  a  single  top  air  gap  of  10 
mils.  A  10  mode  correction  is  presented  in  Figures  5.8  and  5.9. 

The  real  permittivity  extraction  of  Figure  5.8  (a)  represents  an  unsuccessful 
attempt  to  improve  e ',  as  both  corrections  are  less  than  the  truth  data.  This  may 
be  due  to  a  strong  capacitive  effect  at  the  measurement  frequency.  The  corrected 
and  uncorrected  data  sets  corresponding  to  e"  in  Figure  5.8  (b)  are  nearly  identical 
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throughout  the  entire  band,  which  is  consistent  with  previous  measurements.  The 
10  mode  correction,  when  applied  to  magnetic  permeability,  is  very  accurate,  having 
better  performance  than  the  RPI  formulation,  which  can  be  verified  in  Figure  5.9. 
It  must  also  be  pointed  out  that  the  2  mil  uncertainty  in  height  and  thickness  has 
almost  no  effect  on  extraction  of  permeability.  This  is  in  contrast  to  permittivity, 
which  is  clearly  seen  to  be  affected  by  a  small  amount  of  uncertainty  in  the  specified 
sample  dimensions. 

Sample  3,  having  a  top  gap  of  45  mils,  was  also  tested.  A  15  mode  correction 
was  used,  and  the  results  are  displayed  in  Figures  5.10  and  5.11.  The  extraction 
of  e' ,  in  Figure  5.10  (a),  reveals  that  use  of  the  15  mode  correction  gives  a  dramatic 
improvement  over  the  increasingly  worse  uncorrected  data.  However,  even  when  com¬ 
bined  with  the  RPI  formulation,  the  correction  is  not  stable  across  the  entire  band. 
The  standard  mode  correction  would  be  preferred  in  this  case,  since  it  has  a  better 
average  value  and  the  average  RPI  value  is  too  low. 

The  extraction  of  e"  in  Figure  5.10  (b)  is  the  first  instance  when  the  uncorrected 
data  did  not  accurately  approximate  the  truth  data.  In  addition,  the  two  modal 
corrections  show  significant  variation  with  frequency.  This  may  be  due  in  part  to  the 
instability  of  the  Newton-Raphson  Root  search  when  used  in  parameter  extractions 
on  lossless  materials.  It  is  well  known  that  “good”  guesses  are  critical  when  using 
this  algorithm  [1,13].  Of  the  corrections  to  permeability,  shown  in  Figure  5.11,  the 
mode-matching  correction  using  RPI  was  superior  for  both  p!  and  fi" . 


5.2.3  Complex  Propagation  Wave  Number.  It  has  been  mentioned  that  a 
valid  PFW  correction  is  contingent  upon  satisfaction  of  the  transcendental  equation 


k  k 

— sin  ky\h  cos  ky2(b  —  h)  H — —  sin  ky2{b  —  h )  cos  ky\h  =  0 
ei  '  ^2 
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- 10  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


8.5 

8 


7.5 
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8.5  9  9.5  10  10.5  11  11.5  12 

Frequency  (GHz) 

(a) 


- 10  Modes - RPI  - True  (NRW) . Uncorrected 


(b) 

Figure  5.8:  Permittivity  of  FGM-125  in  a  PFW  with  10  mil  top  air  gap.  A  10  mode 
correction  is  given.  Error  bars  indicate  ±2  mil  uncertainty  in  height  and  thickness. 

(a)  Both  corrections  are  unusable,  although  are  preferred  to  the  uncorrected  data. 

(b)  The  modal  correction  and  the  uncorrected  data  are  closest  to  the  truth  data. 
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- 10  Modes  -  -  -RPI  - True  (NRW) . Uncorrected 


- 10  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


Figure  5.9:  Permeability  of  FGM-125  in  a  PFW  with  10  mil  top  air  gap.  A  10 
mode  correction  is  given,  (a)  Non-RPI  modal  correction  is  nearly  exact,  (b)  All  three 
data  sets  are  equally  close  to  the  truth  data. 
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- 15  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


- 15  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


GO 

Figure  5.10:  Permittivity  of  FGM-125  in  a  PFW  with  45  mil  top  air  gap.  A  15 
mode  correction  is  given,  (a)  Although  both  mode  corrections  are  in  the  vicinity  of 
the  truth  data,  the  non-RPI  is  preferred,  (b)  The  RPI  correction  is  the  best  of  the  3 
data  sets,  but  still  varies  significantly  from  the  truth  data. 
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- 15  Modes  -  -  -RPI  - True  (NRW) . Uncorrected 


(a) 


- 15  Modes  -  -  -  RPI  - True  (NRW) . Uncorrected 


Figure  5.11:  Permeability  of  FGM-125  in  a  PFW  with  45  mil  top  air  gap.  A  15 
mode  correction  is  given.  (a,b)  Use  of  RPI  gives  the  best  correction,  as  the  non-RPI 
is  closer  to  the  uncorrected  data. 
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which,  through  finding  7Z,  enables  the  determination  of  the  wave  numbers  ky ±  and 
ky2-  Using  the  procedure  outlined  in  Section  4.2.6,  an  initial  guess  of  7 ffw,  the 
propagation  wave  number  of  a  fully  filled  waveguide,  is  provided  to  the  root  search. 
It  is  trivial  to  calculate  7 ffw  for  any  number  of  desired  modes,  since  it  is  based  on 
the  waveguide  geometry.  However,  using  a  numerical  method  to  determine  7 pfw  may 
lead  to  erroneous  results,  due  to  instability  in  the  iteration. 

Therefore,  it  is  useful  to  observe  the  response  of  the  mode-matching  correction 
as  successive  modes  are  included  in  the  analysis.  This  will  be  done  with  FGM-125 
Sample  3  (S-band),  having  an  air  gap  of  81  miles.  Corrections  to  er  and  /ir  using  1, 
2,  3,  6  and  10  modes  are  illustrated  in  Figures  5.12  and  5.13.  In  Figure  5.12  (a),  it  is 
clearly  seen  that  the  uncorrected  permittivity  obtained  using  NRW  is  far  below  the 
true  value  rendering  it  completely  inaccurate.  Using  the  mode-matching  technique 
with  1  mode  (i.e.  no  higher-order  modes)  yields  a  small  improvement.  Adding  1  higher 
order  mode  improves  the  measurement  further.  Including  2  higher  order  modes  (i.e. 
3  modes  total)  results  in  an  overcorrection,  but  as  more  modes  are  considered,  the 
modal  correction  moves  closer  to  the  true  value.  Similar  observations  can  be  made 
about  the  e"  correction  in  Figure  5.12  (b),  although  in  that  case  the  uncorrected  data 
is  closest  to  the  truth  data.  This  is  also  true  for  the  corrections  to  /T  and  n"  in  Figure 
5.13. 

The  question  remains  as  to  the  reason  for  the  mode-matching  overcorrection 
of  e'  in  Figure  5.12  (a).  Recall  that  the  values  of  the  propagation  wave  number  7 
are  obtained  from  satisfaction  of  the  eigenvalue  equation  (5.5a).  In  order  to  observe 
the  behavior  of  the  eigenvalue  equation  in  the  complex  plane,  two  surface  plots  are 
presented  in  Figure  5.14.  The  material  under  test  is  FGM-125,  observed  at  3.2  GHz, 
using  two  different  air  gap  sizes:  27  mils  and  81  mils,  which  correspond  to  Sample  2 
and  Sample  3  respectively. 

The  ratio  of  a  to  /3,  or  the  propagation  -  attenuation  ratio,  for  each  root  of  7 
from  Figure  5.14  is  given  in  Table  5.4.  A  small  ratio,  such  as  that  for  mode  1  of  the 


69 


6 - 10 - True  (NRW) - Uncorrected 


w 


(a) 


- 1  - 2 - 3  6 - -10  - True  (NRW) - Uncorrected 


Figure  5.12:  Permittivity  of  FGM-125  in  a  PFW  with  81  mil  top  air  gap.  Multiple 
mode  corrections  are  given,  (a)  Too  few  modes  (e.g.  3)  is  an  “overcorrection” ,  but  as 
more  modes  are  included  the  correction  rights  itself,  (b)  Similar  to  (a),  it  is  better  to 
use  more  modes  than  too  few.  Uncorrected  data  gives  the  best  result. 
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6 - 10 - True  (NRW) - Uncorrected 


(a) 

- 1  - 2 - 3  6 - 10  - - True  (NRW) - Uncorrected 


Figure  5.13:  Permeability  of  FGM-125  in  a  PFW  with  81  mil  top  air  gap.  Multiple 
mode  corrections  are  given,  (a)  Minimal  improvement  from  using  increasing  numbers 
of  modes,  but  the  over  correct  ion  described  in  Figure  5.12  occurs  here  as  well,  (b) 
Uncorrected  data  again  gives  the  best  results,  and  is  nearly  identical  to  the  1  &  2 
mode  correction. 
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27  mil  gap,  indicates  the  the  mode  is  highly  propagating,  whereas  a  large  ratio,  such 
as  that  for  mode  5,  indicates  strong  attenuation. 

It  is  expected  that  the  the  dominant  mode  in  a  PFW  system  will  have  a  small 
propagation- attenuation  ratio,  and  this  can  be  verified  in  Figure  5.14  (a),  where  the 
dominant  mode  can  clearly  be  identified.  The  fact  that  the  higher  order  mode  roots 
in  (a)  are  close  to  zero  indicates  minimal  propagation.  This  is  similar  to  the  behavior 
of  modes  in  a  fully  filled  guide,  in  which  the  higher  order  modes  (if  they  even  exist) 
attenuate  rapidly. 

However,  for  the  81  mil  gap  case  in  Figure  5.14  (b),  the  situation  is  very  different. 
Both  the  real  and  imaginary  components  of  the  roots  are  much  larger  than  in  (a).  This 
causes  the  a-(3  ratio  at  every  mode,  except  the  first,  to  be  much  smaller  than  the  27 
mil  gap  case.  The  first  mode  has  7  times  the  attenuation  of  the  27  mil  gap  dominant 
mode,  so  it  can  immediately  be  recognized  that  each  7,  especially  in  the  large  gap 
scenario,  is  both  highly  propagating  and  highly  attenuated.  The  combination  of  these 
effects  implies  that  many  higher  order  modes  are  necessary  for  a  proper  correction 
using  the  PFW  mode-matching  technique.  This  is  consistent  with  the  observations 
made  of  the  successive  mode  corrections  of  Figures  5.12  and  5.13.  Including  several 
higher  order  modes  balances  the  total  held,  which  leads  to  a  more  accurate  correction. 

Table  5.4:  Propagation  -  Attenuation  Ratios 


FGM-125  FGM-125 

27  mil  gap  81  mil  gap 


Mode  # 

a 

P 

a 

p 

a 

P 

a 

p 

1 

11 

40 

0.3 

259 

125 

2.1 

2 

117 

12 

9.8 

325 

109 

3.0 

3 

232 

5 

46.4 

394 

90 

4.4 

4 

347 

3 

115.7 

464 

77 

6.0 

5 

461 

2.8 

164.6 

534 

67 

8.0 

The  potential  exists  for  instability  in  the  numerical  root  search  algorithm  used 
to  converge  to  roots  of  the  eigenvalue  equation.  If  the  solution  space  has  a  shallow 
gradient,  such  in  some  regions  of  Figure  5.14  (b),  the  Newton-Raphson  method  (which 
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uses  a  tangent  line  to  approximate  the  next  guess)  may  obtain  a  poor  approximation. 
It  also  is  helpful  if  the  roots  are  clearly  distinguishable  from  one  another,  as  in  Figure 
5.14  (a).  In  cases  where  the  roots  are  very  close  together,  or  do  not  stand  out  greatly 
from  the  surrounding  space,  it  may  be  appropriate  for  a  different  root  search,  such  as 
the  Muller  Method,  to  be  used  [1,13],  although  this  step  was  not  taken  here. 
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Real  (y) 


(a) 


Figure  5.14:  Two  surface  plots  of  the  eigenvalue  equation  over  the  complex  plane, 
for  which  the  true  PFW  7  will  force  the  output  to  zero.  Material  under  test  is  FGM- 
125,  observed  at  3.2  GHz.  Results  are  presented  in  logarithmic  scale,  (a)  Sample  2, 
27  mil  air  gap,  first  6  PFW  modes.  The  imaginary  part,  or  phase  constant,  is  close  to 
zero,  or  cutoff,  for  the  last  4  modes.  This  ordering  resembles  a  plot  of  7  that  would 
be  expected  in  a  FFW.  (b)  Sample  3,  81  mil  air  gap,  first  5  PFW  modes.  Roots  of  7 
have  large  real  and  imaginary  components,  indicating  both  high  attenuation  and  also 
strong  propagation. 
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VI.  Conclusions  and  Recommendations 


This  thesis  has  demonstrated  the  feasibility  of  using  modal  analysis  to  accurately  per¬ 
form  electromagnetic  material  characterization  in  a  partially  filled  waveguide.  Waveg¬ 
uide  expansion  is  a  common  problem  in  high-temperature  material  measurements, 
since  gaps  form  between  the  sample  and  waveguide  walls.  This  discontinuous  geome¬ 
try  excites  higher  order  which  are  not  accounted  for  in  algorithms  such  as  Nicolson- 
Ross-Weir.  A  correction  was  accomplished  by  considering  a  single  air  gap  between 
the  top  of  the  sample  and  the  waveguide.  Because  the  waveguide  is  partially  filled  in 
the  y- axis,  the  scattered  mode  set  can  be  completely  described  using  TMy  modes. 

Using  Maxwell’s  equations  it  was  shown  that  the  transverse  field  vector  ex  was 
linearly  dependent  on  the  derivatives  of  ey  and  hx.  Therefore,  the  the  boundary 
condition  between  the  empty  and  partially  filled  waveguide  regions  depended  only  on 
explicitly  matching  the  transverse  ey  and  hx  mode  vectors,  although  continuity  of  ex 
can  be  shown  numerically.  A  2-dimensional  Newton-Raphson  root  search  was  used 
to  compute  the  required  value  of  the  complex  propagation  constant  of  each  higher 
order  mode  in  the  partially  filled  waveguide,  which  was  then  used  to  construct  a  finite 
number  of  held  modes  in  the  system.  The  theoretical  scattering  parameters  of  this 
system  are  extracted  and  compared  to  the  experimental  scattering  parameters  in  a 
minimization  equation.  The  top-level  root  search  uses  the  minimization  equations  to 
converge  to  the  true  values  of  relative  permittivity  and  permeability  of  the  sample. 

Measurements  of  acrylic  and  magnetic  shielding  material  in  a  PFW  geometry 
were  performed  at  room  temperature  in  S-band  and  X-band.  For  most  air  gaps, 
acceptable  corrections,  i.e.  within  10%  of  the  true  value,  could  usually  be  achieved 
using  15  modes  or  less.  At  S-band,  modal  corrections  for  air  gaps  greater  than  100 
mils  did  not  converge.  Modal  corrections  on  the  45  mil  gap,  the  largest  gap  at  X-band, 
did  not  converge  using  more  than  10  modes. 

It  was  difficult  to  obtain  useful  modal  corrections  of  the  large  gap  PFW  acrylic 
samples.  Even  when  using  a  single  transmission  measurement  to  extract  permittivity, 
i.e.  a  1-D  search,  the  near-lossless  property  of  the  sample  made  the  root  searching 
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unstable,  and  no  amount  of  modes  used  in  the  correction  would  converge.  Corrections 
to  air  gaps  in  magnetic  shielding  material  measurements  were  usually  more  successful, 
unless  a  significant  gap  (60-100  mils)  was  present.  In  most  cases,  the  algorithm  could 
be  expected  to  compute  the  entire  data  set  (201  points)  in  less  than  20  minutes  using 
a  15  mode  correction. 

The  calculation  of  real  permittivity  (regardless  of  material)  is  more  sensitive  to 
a  top  air  gap  than  either  imaginary  permeability  or  real  and  imaginary  permeability. 
This  is  due  to  the  electric  held  vector  discontinuity  that  exists  across  the  material 
boundary,  which  does  not  exist  in  the  two  other  transverse  held  vectors.  Therefore,  it 
is  reasonable  that  the  mode-matching  technique  will  be  most  effective  at  correcting  the 
parameter  most  corrupted  by  the  presence  of  an  air  gap.  If  the  gap  was  in  the  other 
dimension,  i.e.  left-right,  it  is  expected  that  the  permeability  measurement  would 
suffer.  The  modal  correction  was  usually  able  to  obtain  an  acceptable  correction  for 
e'  when  e"  was  very  small.  The  success  of  modal  corrections  to  nr  depended  on  the 
gap  size  and  frequency  range.  It  was  not  unusual  for  the  “uncorrected”  PFW  data  to 
be  the  most  acceptable  data  set. 

The  addition  of  a  reference  plane  independent  formulation  to  the  mode-matching 
technique  contributed  significant  improvements  to  the  correction.  Foremost,  this  in¬ 
dicates  that  samples  were  not  precisely  aligned  with  the  reference  plane.  Secondly, 
by  using  both  forward  and  reverse  scattering  parameters,  the  effects  of  sample  inho¬ 
mogeneity  on  permittivity  and  permeability  extraction  is  reduced.  Based  on  these 
observations,  it  is  recommended  that  future  material  characterizations  include  a  ref¬ 
erence  plane  independent  solution  method. 

The  experiment  was  not  performed  at  high  temperature,  but  is  a  presented  as 
a  proof  of  concept  technique  for  a  partially  filled  waveguide  environment. 

This  research  presents  a  novel  method  for  characterization  improvements  in 
high  temperature  rectangular  waveguide  measurements.  The  foundation  of  the  single 
air  gap  correction  is  a  mode-matching  technique  of  the  dominant  and  higher  order 
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scattered  TMy  modes.  The  mode-matching  is  then  combined  with  reference  plane 
independence  to  eliminate  possible  phase  shift  errors.  Finally,  the  method  is  applied  to 
measurements  involving  a  magnetic  material.  The  combination  of  these  components 
has  not,  to  the  best  knowledge  of  the  author,  been  presented  before. 

6.1  Future  Work 

It  has  been  observed  that  a  vital  part  of  obtaining  a  correct  solution  using  modal 
analysis  is  the  accurate  determination  of  the  complex  propagation  constant  7 pfw-  If 
an  iterative  search  is  used,  care  must  be  taken  to  ensure  that  “good”  initial  guesses 
are  supplied  to  the  algorithm.  As  an  alternative,  however,  a  root-search  method 
which  can  identify  and  sort  a  desired  number  of  mode  propagation  constants  could  be 
sought.  Use  of  finite  element  analysis,  genetic  algorithms,  or  a  less-sensitive  iterative 
root  search  are  possible  options  for  improvement  in  this  area.  With  more  accurate 
values  of  Jpfw ,  corrections  for  larger  gaps  can  be  calculated. 

The  top  air  gap  PFW  analysis  should  be  combined  with  side  air  gaps  to  better 
simulate  a  high-temperature  situation,  since  it  is  expected  that  all  dimensions  of 
the  waveguide  will  expand.  The  room  temperature  tests  had  a  priori  knowledge  of 
the  dimensions,  specifically  the  height,  of  the  sample  under  test.  This  observation, 
however,  cannot  be  made  at  high-temperature.  It  is  recommended  that  known  thermal 
expansion  coefficients  be  combined  with  the  PFW  analysis  to  obtain  a  reasonably 
accurate  guess  for  the  air  gap  between  the  sample  and  waveguide  while  at  high- 
temperature. 


77 


Appendix  A.  Nicolson-Ross-Weir  Algorithm 

Nicolson,  Ross  [23]  and  Weir  [27]  took  expressions  for  the  SR  reflection  coefficient  and 
SR  transmission  coefficient  and  derived  explicit  formulas  for  the  calculation  of  the 
material  parameters  permittivity  and  permeability  (er,pr).  The  derivation  presented 
here  pertains  to  use  in  a  rectangular  waveguide  where  a  ~  2b. 

Given  the  scattering  parameters  SR  and  SR,  the  reflection  coefficient  R  is  given 
by 

R  =  Q±  VQ 2  -  1  (A.l) 


where  Q  is 


Q  = 


(An)2  —  (A21)2  + 1 


2  s 


(A.2) 


11 


The  choice  of  sign  in  (A.l)  which  forces  |i?|  <  1  is  taken.  Now,  the  transmission 
coefficient  P  is  given  by 


P  = 


A 


21 


1  -  R  ■  An 


(A.3) 


The  impedance  z  can  also  be  found,  as  it  is 


z  = 


1  +  R 
1  -  R 


(A.4) 


In  the  rectangular  waveguide,  the  propagation  constant  7  of  the  dominant,  TEi0 
mode  is 

7o  =  Vlkl  -  kl)  (A. 5) 


Finally,  for  a  material  sample  of  thickness  d,  relative  permittivity  and  perme¬ 
ability  can  be  given  by  the  relations 


/ir 


6r 


In  (P)  ■  z 
7o  d 

(MP,) 


2 


(A.6) 

(A.7) 


78 


Appendix  B.  Integral  Proofs 


The  integrals  of  Chapter  IV  are  shown  here  in  expanded  form. 

Mnp  Aylm  fy) 

J  ’  fydy 

V 


/  mry  pny 

=  /  cos  — —  cos  — — ay 
.  b  b 


/  2n7rVi 
=  /  cos  -^-dy 


b/2,  n  =  p 
0,  n  f  p 


(B.l) 


As  a  special  case,  if  n  =  0  (which  would  represent  the  dominant  mode  in  the  empty 
waveguide),  (B.l)  may  be  simplified  to 

AZ"o v  Ay  10;  fy') 

j“  ^yw  '  fydy 

y 

b 

f  imy,i 

=  /  COS  (0)  COS 


=  J  cos2  (0)  dy 
o 

b,  n  =  p  =  0 
0,  n  f  p 


(B.2) 
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Nnp  Unv  +  'V'r, 


np 


'  np 


Unp  (^l,yln?  ^Py) 


e0cos(ky2  ( b-h ))  cos(k^y)cosP?ydv 
8 !  COS  (ky2h)  J  COS^COS  ft  ^ 


=  2  / cos  ~~  y) y  +  cos  +  y ) ydy 


c 

2 


sin  (feyi  -  2fE)  y  sin  (A;yl  +  2f)  y 


(y  -  ?) 


(y  +  vt  ) 
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2  l  (y  -  f )  (A*  +  f ) 


(B.3) 
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